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ABSTRACT 


A/u .i-hr  r 

-  ;This  study  brings  together  the  programs  developed  in  three  previous 
reports  into  a  unified  sampling  system.  One  segment  contains  tne  single  rsno 

I 

double  sampling  plans  of  MIL-STD-105D  along  with  subroutines  that  evaluate  the 
normal,  tightened  and  reduced  pl^ns  including  switching  rules.  Output 
includes  system  OC  curves,  and  A5N,  AOQ,  and  -Afi-  curves  in  either  graphical  or 
tabular  form.  Another  segment  derives  single  and  double  sampling  plans  to 
satisfy  two  points  designated  on  the  OC  curve.  An  algorithm  is  employed  which 
seeks  to  minimize  the  ASN  at  a  designated  Acceptable  Quality  Level  (AOL).  The 
third  segment  is  designed  for  use  when  rectifying  inspection  is  employed.  It 
derives  single  and  double  sampling  plans  based  on  one  point  on  the  OC  curve. 

either  the  process  average  or  a  specified  AOL  value,  intended  to  not  exceed  a 

7.  -■‘1  ?  -.7  r 

specified  AOQtr  and  minimizing  the  AFI  at  the  process  average  or  specified 
AQL.  Programming  is  written  in  FORTRAN  IV  and  development  was  on  a  VAX  11-750 
computer. 


INTRODUCTION 


This  report  brings  together  the  results  of  three  previous  studies  dealing 
with  acceptance  sampling  using  attributes  data.  The  previous  studies  dealt 
with  single  and  double  acceptance  sample  plans.  The  first  involved  a 
computerized  version  of  MIL-STD-105D.  The  second  dealt  with  the  design  of 
double  sampling  plans  in  cases  where  rectifying  inspection  is  not  employed. 

An  objective  function  is  introduced  which  minimizes  the  average  amount  of 
inspection  when  the  process  is  operating  at  what  is  considered  an  acceptable 
level.  The  third  study  assumes  that  rectifying  inspection  is  planned  based  on 
some  specified  Average  Outgoiing  Quality  Limit  and  has  as  an  objective  function 
the  minimization  of  inspection  when  the  process  is  operating  at  an  acceptable 
level.  All  computer  programming  is  in  ANSI  FORTRAN  IV  and  the  programs  were 
developed  on  a  VAX  11-750  computer  utilizing  the  VMS  operating  systems. 

DESCRIPTION  OF  THE  PROGRAM  SEGMENTS 
List  of  Program  Segments 

The  programming  system  is  composed  of  seven  main  segments.  Each  performs 
a  specific  activity  with  respect  to  the  operation  of  the  system  or  contains 
the  code  for  executing  the  user's  choice  of  sampling  plans.  These  include: 

1.  Finding  and  evaluating  MIL-STD-105D  single  and  double  sampling  plans 
as  a  system  of  normal,  tightened,  and  reduced  inspection  plans. 

2.  Deriving  sampling  plans  to  minimize  the  Average  Sample  Number  (ASN), 
or  average  sample  size; 

3.  Deriving  sampling  plans  to  minimize  the  Average  Fraction  Inspected 
(AFI)  based  on  a  specified  Average  Outgoing  Quality  Limit  ( AOQL) . 


The  program  segments  are: 


QUALITY. COM  A  command  file  that  complies  and  links  the  program 

files  into  an  executable  file,  QMAIN.EXE. 

QMAIN.FOR  The  main  program  control  file  from  which  are  called  the 
three  sampling  plan  programs. 

QMIL.FOR  The  program  containing  MIL-STD-105D  sampling  plans  and 
the  analysis  package. 

QASN.FOR  The  program  for  finding  single  and  double  sampling 
plans  based  on  minimizing  the  ASN. 

QAFI.FOR  The  program  for  finding  single  and  double  rectifying 
inspection  sampling  plans  based  on  minimizing  the  AFI. 

PROBS1.FOR  A  program  for  calculating  probabilities,  binomial  or 
Poisson,  for  single  sampling  plans.  It  is  used  by  QASN.FOR  and 
QAFI.FOR. 

PR0BD1.F0R  A  program  for  calculating  probabilities,  binomial  or 
Poisson,  for  double  sampling  plans.  It  is  used  by  QASN.FOR  and 
QAFI.FOR. 


Program  Segment  QMIL.FOR 

The  program  QMIL.FOR  contains  all  of  the  single  and  double  sampling 
acceptance  plans  from  MIL-STD-105D  for  normal,  tightened,  and  reduced 
inspection.  They  are  organized  such  that  the  user  goes  through  identically 
the  same  procedure  as  if  the  Standard  were  being  used.  First,  an  Inspection 
Level  is  entered  (Special  Levels  S-l  through  S-4  or  General  Levels  I,  II,  or 
III)  followed  by  lot  or  batch  size.  This  leads  the  program  to  the  selection 
of  the  appropriate  sample  size  code  letter.  The  required  AQL  (Acceptable 
Quality  Level)  is  entered  next.  This  value  must  be  one  of  the  prescribed 
values  in  the  Standard.  Any  other  value  will  lead  to  extraneous  and  undesired 
results.  Those  values  are,  expressed  in  percent  defective  or  defects  per  100 
units: 

0.010,  0.015,  0.025,  0.040,  0.065,  0.10,  0.15,  0.25,  0.40,  0.65, 

1.0,  1.5,  2.5,  4.0,  6.5,  10,  25,  40,  65,  100,  150,  250,  400,  650,  1000. 
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Once  a  system  of  sampling  plans  has  been  obtained,  the  user  is  qiven  the 
option  of  having  the  system  analyzed.  Output  of  the  analysis  may  be  either  in 
tabular  or  graphical  form.  In  tabular  form,  a  range  of  values  of  p  for  which 
the  calcuations  are  to  be  made  must  be  entered.  If  the  user  has  no  idea  of 
the  range  of  p  required,  it  may  be  desirable  to  first  obtain  graphical  output 
and  then  select  the  values  of  p  for  which  accurate  results  are  to  be 
obtained.  One  caution  with  respect  to  graphical  output  is  necessary  at  this 
point.  In  order  to  obtain  clear,  easily  readable  graphical  output,  the 
printer  unit  must  be  capable  of  and  set  on  condensed  printing.  Graphical 
output  uses  more  space  than  the  standard  80  column  pica  type  font  allows. 
Hence,  when  restricted  to  80  columns,  two  lines  are  used  for  each  actual  line 
of  output. 

Output  of  the  analysis  in  tabular  form  appears  as  follows: 

1.  The  value  of  p  for  which  the  calculation  was  made. 

2.  The  probability  of  acceptance. 

3.  The  average  sample  size  (ASN). 

4.  The  average  outgoing  quality  (AOQ). 

5.  The  average  fraction  inspected  (AFI). 

Item  3  applies  only  to  those  items  in  samples.  Items  4  and  5  apply  only  if 
rectifying  inspection  is  being  employed;  i.e.,  that  all  rejected  lots  are 
screened  of  defective  units  (defects)  and  are  repaired/ replaced  in  the  lot  in 
acceptable  condition. 

The  graphical  output  provides  the  same  information  except  that  curves  of 
the  characteristics  are  plotted  rather  than  numerical  values  printed.  Because 
of  limitations  on  the  type  of  graphical  output,  the  tabular  output  will  be 
more  accurate. 

Analytical  results  of  the  program  will  not  correspond  exactly  with  those 
shown  in  MIL-STD-105D  because  the  program  treats  the  normal,  tightened,  and 
reduced  inspection  plans,  along  with  the  switching  rules  for  changing  from  one 


to  the  other,  as  a  system  of  plans  and  analyzes  the  entire  system  rather  than 
each  plan  individually.  The  system  is  analyzed  as  a  Markov  chain.  Steady 
state  probabilities  of  being  in  each  of  a  series  of  states  under  normal, 
tightened  and  reduced  inspection  are  calculated  for  each  value  of  p  in  the 
effective  range  of  the  OC  curve.  These  probabilities  are  then  combined  to 
give  probabilities  of  being  in  normal,  tightened  and  reduced  inspection.  From 
these  values  are  calculated  the  probabilities  of  acceptance  of  lots  or 
batches,  the  OC  curve  points,  the  ASN,  AFI,  and  AOQ.  As  is  the  case  in  the 

Standard  itself,  each  calculation  assumes  that  lots  are  being  formed  from  a 

process  generating  a  constant  value  of  p,  i.e.,  a  process  in  statistical 
control . 

Details  of  the  models  and  procedures  are  contained  in  reference  (1). 

Program  Segment  QASN.FOR 

This  program  segment  enables  the  user  to  design  custom  double  sampling 
plans  to  minimize  the  average  sample  size  based  on  two  points  on  the  operating 

characteristic  curve.  In  the  process,  it  also  specifies  the  minimum  single 

sampling  plan  that  meets  the  OC  curve  requirements.  The  selection  of  plans 
based  on  ASN  minimization  are  intended  for  use  when  rectifying  inspection 
either  is  not  or  can  not  be  used.  Inputs  to  this  program  segment  include  a 
choice  between  the  binomial  and  Poisson  distributions  for  calculating  proba¬ 
bilities  of  acceptance,  the  two  specified  points  on  the  OC  curve,  a  seed  value 
for  the  sample  size  of  the  equivalent  single  sampling  plan,  and  a  value  for 
the  rejection  number  on  the  first  sample  of  the  double  sampling  plan. 

The  output  consists  of  the  minimum  sampling  plan,  a  series  of  double 
sampling  plans,  the  double  sampling  plan  found  to  have  the  minimum  ASN 
evaluated  at  the  input  good  quality  level,  and  the  maximum  value  for  the  ASN 
for  the  R1  value  orginally  input  to  the  program. 


The  two  points  designated  for  the  OC  curve  are  related  to  the  traditional 


Producer's  and  Consumer's  risks.  The  program  uses  PO  to  designate  the 
producer's  maximum  acceptable  quality  level  (an  AQL  as  defined  in  MIL-STD- 
105D)  or  process  average  with  the  minimum  probability  of  acceptance  set  at 
(1-  «) .  In  this  case,  a  is  the  designated  maximum  producer's  risk.  The 
quality  level  PI  is  the  designated  consumer's  maximum  acceptable  quality  level 
(RQL,  or  rejectable  quality  level)  with  a  maximum  probability  of  acceptance 
of  8.  In  this  case,  8  is  the  designed  maximum  Consumer's  risk.  The  equations 
are  of  the  form: 

L(P0)  >  1  -  a 
L  ( PI )  <  8 

where  L(P)  is  the  likelihood  function  (probability  of  acceptance  formula)  for 
the  particular  sampling  plan  being  analyzed. 

With  these  two  equations  alone,  an  infinite  number  of  sampling  plans, 
single  and  double,  exist  which  satisfy  the  inequality  constraints.  A 
reasonable  choice  from  among  that  group  is  that  plan  which  provides  for 
minimum  inspection  on  the  average.  Thus  an  objective  function  is  introduced 
that  minimizes  the  ASN  when  the  process  is  operating  at  or  below  the  AQL. 

This  value  is  designated  ASN  or  ASN(PO)  on  the  computer  output. 

Users  familiar  with  MIL-STD-1050  will  remember  that  double  and  multiple 
sampling  plans  in  that  standard  are  designed  such  that  their  respective  ASN's 
never  exceed  the  sample  size  of  the  single  sampling  plan  with  the  same  OC 
curve.  The  capability  of  establishing  this  additional  constraint  has  been 
provided  for  in  this  program  segment.  As  part  of  the  output,  the  user 
receives  the  maximum  value  of  the  ASN  reached  for  each  plan  analyzed.  This 
values  is  designed  ASNMAX.  (The  program  searches  the  ASN  function  as  a 
function  of  p  and  outputs  that  value.)  Thus  the  user  may  search  through  the 
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program  output  for  those  plans  with  a  ASNMAX  less  than  or  equal  to  the  minimum 
single  sampling  plan  and  select  from  them  the  plan  with  the  minimum  ASN(PO). 

It  is  worth  noting  at  this  point  that,  in  most  of  the  plans  explored  in  the 
development  of  this  segment,  addition  of  this  constraint  did  not  result  in 
changing  the  optimal  choice. 

Reference  (2)  presents  detailed  findings  and  a  full  description  of  the 
program's  development. 

Program  Segment  QAFI.FOR 

This  program  segment  develops  single  and  double  sampling  plans  which  meet 
a  specified  Average  Outgoing  Quality  Limit  (AOQL).  They  are  intended  for  use 
only  when  rectifying  inspection  (100%  inspection  of  the  balance  of  rejected 
lots)  is  employed.  Calculations  assumes  that  rectifying  inspection  is  100% 
effective  and  that  the  binomial  distribution  provides  sufficient  accuracy  in 
calculating  probabilities  of  acceptance. 

Inputs  to  this  segment  include  a  Producer's  Risk  Point  (PO,  1-a),  the 
desired  AOQL  (designated  AOQL*),  and  the  lot  or  batch  size.  Outputs  include 
the  minimum  single  sampling  plan  satisfying  the  constraints,  a  series  of 
double  sampling  plans  all  of  which  satisfy  the  AOQL  constraint,  and  that  plan 
which  satisfies  the  constraints  and  minimizes  the  Average  Fraction  Inspected 
(AFI)  when  the  process  is  operating  at  the  quality  level  PO.  The  quality 
level  PO  may  be  interpreted  either  as  an  AQl  value  as  defined  in  MIL-STD-105D 
or  as  a  process  average  as  defined  and  employed  in  the  Dodge-Romig  tables.  As 
in  the  Dodge-Romig  tables,  the  rejection  numbers  on  the  first  and  second 
samples,  R1  and  R2,  are  set  at  the  second  acceptance  number,  C2,  plus  one. 
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The  objective  function  and  constraint  equations  are  of  the  form: 


minimize  AFI(PO) 

subject  to: 

L(P0)  >  1  -  a 
AOQL  <  AOQL* 

R1  =  R2  =  C2  +  1 

First,  the  program  finds  the  minimum  single  sampling  plan  satisfying  the 
likelihood  function  and  AOQL  constraints.  This  information  is  used  internally 
to  assist  in  setting  bounds  for  the  algorithm  that  solves  for  the  double 
sampling  plans.  Then  the  algorithm  shifts  to  seeking  the  double  sampling  plan 
satisfying  the  constraints.  In  its  search,  a  number  of  plans  are  found  and  it 
is  from  this  group  that  the  plan  that  minimizes  AFI(PO)  is  selected.  Output 
includes  the  sample  sizes  Nl,  and  N2,  acceptance  numbers  Cl  and  C2,  and  the 
AFI(PO),  designated  simply  as  AFI. 

Reference  (3)  provides  a  complete  description  and  analysis  of  the 
development  and  operation  of  the  algorithm  and  program. 

PROGRAM  OPERATING  INSTRUCTIONS  AND  OUTPUT 
Initializing  the  Program 

In  the  following  instructions,  user  inputs  are  in  capital  letters.  A 
backward  arrow  (-*-)  indicates  a  carriage  return. 

Program  initialization  begins  by  entering: 

@QUALITY<- 

This  instruction  executes  the  command  file  (OUALITY.COM)  which  compiles  and 
links  the  program  files  into  a  single  executable  file  (QMAIN.EXE).  If 
QMAIN.EXE  has  been  created  and  exists  in  the  file,  this  step  may  be 
el iminated. 

Once  QMAIN.EXE  has  been  created,  the  user  enters: 


RUN  QMAIN+ 

This  command  starts  program  execution.  Very  shortly  the  following  message 
will  appear  on  the  terminal  screen: 


WHAT  DO  YOU  WISH  TO  DO? 

1- DERIVE  SAMPLING  PLANS  TO  MINIMIZE  ASN 

2- DERIVE  SAMPLING  PLANS  TO  MINIMIZE  AFI 

3- EVALUATE  MIL-STD-105D  SAMPLING  SCHEME 

4- EXIT  THIS  PROGRAM 

Entering  the  number  for  your  selection  followed  by  a  carriage  return  transfers 
program  control  to  QASN.FOR,  OAFI.FOR,  QMIL.FOR,  or  the  computer  operating 
system.  These  instructions  are  repeated  after  each  run  of  the  program  until 
the  user  exits  with  instruction  (4). 


Running  QASN.FOR 

Entering  number  1  from  the  main  transfers  program  control  to  QASN.FOR, 

the  program  for  generating  double  sampling  plans  designed  to  minimize  the  ASN 

at  quality  level  PO.  What  follows  is  the  menu  QASN.FOR. 

The  first  question  in  this  menu  asks  you  to  name  an  output  file  (8 

alphanumeric  characters  or  less). 

WHAT  DO  YOU  WISH  TO  CALL  YOUR  OUTPUT  FILE? 

User  enters  a  name  and 

Subsequent  questions  are: 

CODES  FOR  SELECTING  APPR.  PROB.  DIST. 

BINOMIAL  =1 

POISSON  =2 

User  then  enters  the  appropriate  code.  The  Poisson  uses  the  intensity 
parameter  A  =  np  in  terms  of  defects  per  100  units. 

SELECT 

SAMPLE  PLANS  ONLY  =  1 
ASN  VALUES  ONLY  =  2 
OR  BOTH  =  3 


8 


The  most  useful  of  these  choices  is  to  enter  3.  Entering  1  produces  an  output 
of  a  large  number  of  sampling  plans  all  of  which  meet  th  criteria  for  the  two 
points  on  the  OC  curve.  Output  includes  the  values  of  Cl,  C2,  and  Nl,  the 
range  over  which  N2  may  operate,  and  the  probabilities  of  acceptance  at  PI  and 
PO,  respectively.  Entering  2  results  in  the  production  of  a  large  number  of 
plans  satisfying  the  OC  curve  constraints  and  lists  values  of  ASNMAX  and 
ASN(PO).  When  3  is  entered,  the  optimization  algorithm  is  called  in  thus  the 
output  is  limited  to  those  plans  which  were  considered  candidates  for  the 
optimum.  Output  includes  Cl,  C2,  and  Nl,  the  effective  range  for  N2,  and  the 
ASNMAX  and  ASN(PO)  when  the  smallest  value  of  N2  is  used.  It  then  selects  the 
global  optimum  plan  for  the  input  value  of  Rl. 

INPUT  ALPHA  User  enters  the  Producer's  Risk,  <*-*- 

INPUT  BETA  User  enters  the  Consumer's  Risk,  B<- 

INPUT  PO  User  enters  P0<- 

INPUT  PI  User  enters  Pl<- 

INPUT  A  SEEO  FOR  THE  SINGLE  SAMPLING  NO. 

IF  NO  SEEO  AVAILABLE  ENTER  ZERO  AS  THE  SEED  VALUE 

If  the  user  has  some  idea  what  the  sample  size  will  be  for  the  single  sampling 
plan  satisfying  the  OC  curve  requirements,  a  conservative  value  somewhat  lower 
than  that  number  may  be  entered.  This  option  is  useful  only  when  the  sample 
size  is  known  to  be  large  and  user  wishes  to  save  some  time.  Its  use  may  lead 
to  nonoptimal  results  if  the  seed  value  is  larger  than  the  optimal  value. 

INPUT  A  VALUE  FOR  (Rl  -  C2) 

IF  Rl  =  C2  THEN  THE  VALUE  WOULD  BE  0 

IF  Rl  >  C2  THEN  THE  VALUE  WOULD  BE  A  POSITIVE  NO. 

IF  Rl  <  C2  THEN  THE  VALUE  WOULD  BE  A  NEGATIVE  NO. 

This  option  allows  the  user  to  vary  the  value  of  Rl  in  the  search  for 
optimality.  Each  run  of  the  program  provides  the  optimal  plan  for  a  given 


selected  value  of  Rl. 


From  this  point  on  the  program  takes  over  generating  the  desired  output. 


Running  QAFI.FOR 

Entering  2  from  the  main  menu  transfers  control  to  QAFI.FOR,  the  program 
segment  for  finding  rectifying  inspection  plans  satisfying  AOQL  and  AQL 
constraints  while  minimizing  the  AFI  at  the  AQL  level.  The  menu  for  QAFI.FOR 
leads  the  user  as  follows: 

ENTER  VALUE  OF  ALPHA  User  enters  Producer's  Risk,  CH¬ 
ESTER  VALUE  OF  PO  User  enters  quality  level  PO*- 

ENTER  AOQL  VALUE  User  enters  AOQL*-  (in  percent) 

ENTER  LOT  SIZE  User  enters  lot  size*- 

The  program  then  proceeds  to  generate  a  series  of  sampling  plans  for 
which  R1  =  R2  =  C2  +  1  and  selects  the  optimal  plan  from  among  this  group. 
Output  includes  the  values  of  Cl,  C2,  Nl,  N2,  and  AFI(PO),  respectively,  and 
the  optimal  plan. 

Running  QMIL.FOR 

Entering  3  from  the  main  menu  transfers  control  to  QMIL.FOR,  the  program 

segment  for  finding  single  and  double  sampling  plans  from  MIL-STD-105D  and,  at 

the  user's  discretion,  evaluating  the  resulting  normal,  tightened  and  reduced 

inspection  plans  as  a  system.  Output  of  the  evaluation  may  be  either 

graphical  or  tabular.  The  menu  for  this  segment  is  as  follows: 

WHAT  00  YOU  WISH  TO  CALL  YOUR  OUTPUT  FILE? 

User  enters  a  name  for  the  output  file*-  (8  alphanumeric  characters  or  less). 

ENTER  INSPECTION  LEVEL  IN  QUOTATION  MARKS 
E.G.,;  SPECIAL:  'SI' ,  'S2'  ,  'S3'  ,  'S4‘ 

GENERAL:  '1',  '2',  '3' 

User  enters  the  appropriate  code  followed  by  *-.  A  note  of  caution  here.  Any 
error  on  entry  will  cause  program  failure  and  automatic  exit.  Depending  on 
the  terminal  being  used,  it  may  be  necessary  to  use  apostrophe  marks  rather 
than  quotation  marks  surrounding  the  entry. 


.  - 
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ENTER  LOT  SIZE:  User  enters  lot.  size'*- 

ENTER  AQL  IN  PERCENT.  REMEMBER.  ONLY  A  STANDARD 
AQL  IS  ALLOWABLE  User  enters  AOL  «■ 

DO  YOU  WANT  SINGLE  (‘S')  OR  DOUBLE  ( ‘ D* ) 

SAMPLING  PLANS;  (NOTE:  ENTER  S  OR  D  IN  QUOTES) . 

The  same  caution  for  entry  of  the  inspection  level  applies  here  as  well. 

The  program  then  proceeds  to  select  the  desired  sampling  system  from  the 

Standard.  Once  the  plans  are  displayed,  user  is  asked: 

00  YOU  WANT  SCHEME  EVALUATION..? 


IF  YES  ENTER . 1 

IF  NO  ENTER . 2 


If  the  answer  to  this  question  is  no,  program  execution  stops  and  the  user  is 
presented  with  the  main  menu.  If  the  answer  is  yes,  execution  continues  with: 
DO  YOU  WANT  A  TABLE  OR  A  GRAPH  FORMAT? 

FOR  GRAPH . ENTER:  1 

FOR  TABLE . ENTER:  2 

Entering  l-*-  for  graphical  format  results  in  the  output  of  a  plot  of  the  OC 
curve,  the  ASN  curve,  the  AOQ  curve,  and  the  AFI  curve.  (Naturally,  the  A00 
and  AFI  curves  have  no  meaning  unless  rectifying  inspection  is  intended.)  In 
order  for  these  curves  to  be  readable,  the  printer  must  be  capable  of  and  set 
on  compressed  printing.  The  horizontal  scale  of  each  of  the  graphs 
(probability  scale  in  the  case  of  the  OC  curve)  requires  more  than  80  columns. 
Thus,  if  an  attempt  is  made  to  print  graphs  in  standard  80  column  format,  the 
printing  will  occur  on  two  consecutive  lines  making  the  whole  thing  look  very 
weird. 

When  2-*-  is  entered,  i.e.,  a  tabular  format  is  requested,  the  computer  returns 
with  the  following  questions: 
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SPECIFY  THE  NUMBER  OF  FRACTION  DEFECTIVE  VALUES.. 

The  user  must  enter  the  number  of  values  of  p  for  which  the  equation  is  to  be 
made. 

ENTER  THE  FRACTION  DEFECTIVE  VALUE(S), 

(PUT  A  COMMA  BETWEEN  VALUES.) . 

The  user  then  enters  the  values  as  fractional  quatntities,  not  percents.  In 
this  case  the  output  is  a  table  containing  p,  P(A),  ASN,  AOQ  and  AFI, 
respectively.  It  may  be  useful  to  obtain  a  graphical  output  on  a  first  run  to 
obtain  the  effective  range  of  p  required  for  a  good  tabular  output. 

Naturally,  the  table  is  more  accurate  than  is  the  graph. 
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APPENDIXES 


Appendix  A  contains  several  example  runs  for  each  of  the  program  segments. 
The  entire  program  is  listed  in  Appendix  B. 


Example  1  QASN.FOR 


WHAT  DO  YOU  WISH  TO  DO? 

1 - DERIVE  PLANS  TO  MINIMIZE  ASN 

2—  DERIVE  PLANS  TO  MINIMIZE  AFI 

3—  EVALUATE  MIL-STD-105D  SAMPLING  SCHEME 

4- EXIT  THIS  PROGRAM 


ENTER  CHOICE 
-i- 

WHAT  DO  YOU  WISH  TO  CALL  YOUR  OUTPUT  FILE? 
PR1 


CODES  FOR  SELECTING  APPR.  PROB.  DIST. 


BINOMIAL 

POISSON 


-1 

•2 


SELECT 


SAMPLE  PLANS  ONLY  -I 
ASN  VALUES  ONLY  -2 
OR  BOTH  -3 


.05 

INPUT 

ALPHA 

INPUT 

BETA 

.  015 

INPUT 

P0 

.05 

INPUT 

PI 

U 

INPUT  A  SEED  FOR 

0 


IF  NO  SEED  AVAILABLE  ENTER  ZERO  AS  THE  SEED  VALUE 
INPUT  A  VALUE  FOR  (R1-C2) 

IF  R1-C2  THEN  THE  VALUE  WOULD  BE  0 

IF  R1 >  C2  THEN  THE  VALUE  WOULD  BE  A  POSITIVE  NO. 

IF  Rl <C2  THEN  THE  VALUE  WOULD  BE  A  NEGATIVE  NO. 
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1 


ALPHA  -0.0500  BETA  -0.1000 

P0  -0.0150  PI  -0.0500 

REJECTION  NO.  OF  FIRST  SAMPLE  (Rl)  -  C2+ (  1) 


♦SINGLE  SAMPLING  PLAN* 
ACCEPTANCE  NO.  <C>  -  6 
LOWER  BOUND  ON  N  (NS) 
UPPER  BOUND  ON  N  <NL> 

-  209 

-  220 

DOUBLE 

SAMPLING 

PLANS 

FOR  Cl*  0  C2* 

6 

(Nl) 

(N2S) 

(N2L) 

ASNMAX 

ASN 

59 

60 

166 

164 

166 

165 

212. 4477 
211.5781 

156.9425 

157.7704 

DOUBLE 

SAMPLING 

PLANS 

FOR  Cl 

*  1  C2* 

6 

(Nl) 

(N2S) 

(N2L ) 

ASNMAX 

ASN 

89 

90 

139 

136 

139 

137 

201. 8087 
200. 3593 

142. 6555 
143. 2139 

DOUBLE 

SAMPLING 

PLANS 

FOR  Cl 

-  2  C2- 

6 

(Nl) 

(N2S) 

(N2L) 

ASNMAX 

ASN 

115 

116 

US 

112 

116 

114 

191. 8224 
190. 8090 

143. 4053 
144. 1031 

DOUBLE  SAMPLING  PLAN8 


FOR 

Cl-  3  C2- 

6 

<N1 ) 

(N2S) 

(N2L) 

139 

94 

96 

140 

91 

94 

DOUBLE  SAMPLING 

PLANS 

FOR 

Cl-  4  C2- 

6 

(Nl) 

<N2S) 

(N2L) 

161 

82 

82 

162 

76 

80 

DOUBLE  SAMPLING 

PLANS 

FOR 

Cl-  5  C2- 

6 

(Nl) 

(N2S) 

(N2L) 

184 

75 

87 

185 

57 

80 

♦SINGLE  SAMPLING  PLAN# 
ACCEPTANCE  NO.  <C)  =  7 
LOWER  BOUND  ON  N  (NS) 
UPPER  BOUND  ON  N  (NL) 

DOUBLE  SAMPLING  PLANS 


FOR  Cl-  0  C2—  7 


(Nl) 

(N2S) 

(N2L) 

50 

226 

228 

51 

220 

226 

52 

215 

225 

53 

211 

223 

ASNMAX  ASN 

186.4893  153.2870 

185.9685  154.0707 


ASNMAX  ASN 

188.3876  167.9648 

187.3810  168.5728 


ASNMAX  ASN 

196.2480  186.8779 

194.3076  187.2266 


/ 


»  234 
-  266 


ASNMAX  ASN 

265.0669  169.8499 

260.3219  169.2180 

256.5318  169.0233 

253.6964  169.2884 
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DOUBLE  SAMPLING  PLANS 


FOR  Cl-  1  C£«  7 


(N1  ) 

<N2S> 

<  N2L ) 

ASNMAX 

ASN 

SI 

198 

202 

252. 9615 

148.9811 

82 

192 

200 

248.  7259 

148. 9649 

S3 

187 

198 

245. 3612 

149.2348 

DOUBLE  SAMPLING 

PLANS 

FOR 

Cl-  2  C2- 

7 

<N1 ) 

<N2S) 

(N2L) 

ASNMAX 

ASN 

108 

178 

181 

242.  3765 

147. 3077 

109 

170 

179 

237. 3196 

147. 2064 

110 

164 

177 

233.  7746 

147.5018 

DOUBLE  SAMPLING 

PLANS 

FOR 

Cl-  3  C2- 

7 

(Nl) 

(N2S) 

(N2L) 

ASNMAX 

ASN 

134 

157 

164 

231.0574 

156. 3438 

135 

147 

162 

225.  8645 

156. 3163 

136 

139 

160 

221. 9088 

156. 5329 

DOUBLE  SAMPLING 

PLANS 

FOR 

Cl-  4  C2- 

7 

<N1> 

<N2S) 

(N2L  > 

ASNMAX 

ASN 

159 

141 

156 

224.  8282 

171.6242 

160 

125 

152 

218. 3522 

171.4111 

161 

115 

149 

214.6785 

171.7019 

DOUBLE  SAMPLING  PLANS 


FOR  Cl-  5  C2-  7 


(N 1  ) 

(N2S) 

(N2L  > 

ASNMAX 

ASN 

184 

117 

181 

220. 1528 

190. 2280 

183 

96 

172 

214. 6608 

190. 2104 

186 

83 

165 

212. 2600 

190. 7030 

DOUBLE 

SAMPLING 

PLANS 

FOR  Cl 

-  6  C2  = 

7 

( N 1  ) 

(N2S) 

<N2L ) 

ASNMAX 

ASN 

209 

73 

266 

220. 3671 

210. 8894 

210 

57 

266 

218. 6382 

21 1. 4632 

GLOBAL 

MINIMUM 

ASN(P0> 

-  142.66 

CORRESPONDING  N1  «  89 

CORRESPONDING  N2S  -  139 

CORRESPONDING  Cl  -  1 

CORRESPONDING  C2  -  6 
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Example  2  QAFI.FOR 


UHH  iK)  vr  ■  ■  WISH  :0  00” 

DERIVE  Pt  AnF  TD  MINIMIZE  ASN 

_  UfcRIVE  Plans  to  minimize  afi 

3- EVALUATE  MIL-STD-105D  SAMPLING  SCHEME 

4- EXIT  THIS  PROGRAM 


ENTER  CHOICE 
£ 

ENTER  VALUE  OF  ALPHA 


.  05 


EiviER  vhi_UE  OF 


.  015 
.  0g 
£000 


I 


ENTER  AOQL  VALUE 

ENTER  LOT  SIZE 

ALPHA=  0. 050000 
P0=  0.015000 

AOQL=  0. 020000 
N=  2000 


***  SINGLE  SAMPLING  PLAN  *** 


I 


NS=1£0 
C=  4 
AFI (P0) 


«=0.  093093 


**  DOUBLE  SAMPLING  PLANS  ** 


Cl 

C2 

N 1 

N2 

AFI 

0 

4 

£9 

91 

0. 056349 

1 

4 

55 

6G 

0. 060894 

Cl 

C£ 

N 1 

N2 

AFI 

0 

1 


5 

5 


28 

55 


120 

93 


0. 051799 
0. 054159 


0 

1 


6 

6 


147 

12® 


0. 050792 
0. 050957 


Cl 

C2 

N 1 

N2 

AFI 

0 

7 

28 

1  74 

0.051735  { 

1 

7 

55 

147 

0.  049855 

CL 

7 

82 

121 

0. 056838  < 

Cl 

C£ 

N  1 

N2 

AFI 

0 

8 

28 

201 

0.  053922 

1 

8 

55 

174 

0.050058  1 

2 

8 

82 

148 

0. 055690 

SAMPLING 

PLAN 

MIN I MUMS 

Cl=  1  C2=  7 

N 1 =  55  N2= 1 47 

MINIMUM  AFI=0. 049855 


Example  3  QMIL.FOR 


WHAT  DO  YOU  WISH  TO  DO? 

1- DERIVE  PLANS  TO  MINIMI  ZE  ASN 

2- DERIVE  PLANS  TO  MINIMIZE  AF I 

3- EVALUATE  MIL-STD-105D  SAMPLING  SCHEME 

4- EXIT  THIS  PROGRAM 


ENTER  CHOICE 
_3_ 

ENTER  INSPECTION  LEVEL  IN  QUOTATION  MARKS 
E.G.  ;  SPECIAL  :  ’  S 1  *  ,  ’  S3’  ,  ’  S3’  ,  1  S4’ 

GENERAL  :’  1’  ,  ’  2’  ,  •  3’ 

1  2’ 

ENTER  LOT  SIZE: 

2000 

ENTER  AQL  IN  PERCENT.  REMEMBER,  ONLY  A 

STANDARD  AQL  IS  ALLOWABLE 

1.5 

DO  YOU  WANT  SINGLE (’S’)  OR  DOUBLE  < ’ D*  > 
SAMPLING  PLANS  ; (NOTE : ENTER  S  OR  D  IN  QUOTES). 
’  S’ 


THESE  PLANS  ARE: 


*: NORMAL: ::::::::::: : TIGHTENED: :::::::::: : REDUCED: :* 


AC 

1=  5 

AC 

1  = 

3 

AC  1-  2 

RE 

1=  6 

RE 

1  = 

4 

RE  1-  5 

N= 

125 

N= 

125 

N=  50 

DO  YOU  WANT  SCHEME  EVALUATION..? 

IF  YES  ENTER . 1 

IF  NO  ENTER . 2 

2 


Example  4  QMIL.FOR 

WHAT  DO  YOU  WISH  TO  DO? 

1 - DERIVE  PLANS  TO  MINIMIZE  ASN 

2-  DERIVE  PLANS  TO  MINIMIZE  AFI 

3- EVALUATE  MIL-STD-105D  SAMPLING  SCHEME 

4- EXIT  THIS  PROGRAM 


ENTER  CHOICE 
3 

ENTER  INSPECTION  LEVEL  IN  QUOTATION  MARKS 
E.G.  ;  SPECIAL  : 1  SI  ’  ,  ’  S2’  ,  ’  S3’  ,  •  S4’ 

GENERAL  :»1»  ,  ’  2’ , ’ 3’ 

LSI 

ENTER  LOT  SIZE: 

2000 

ENTER  AQL  IN  PERCENT.  REMEMBER,  ONLY  A 

STANDARD  AQL  IS  ALLOWABLE 

1.5 

DO  VOU  WANT  SINGLE  <  ’  S’  )  OR  DOUBLE(*D’> 
SAMPLING  PLANS  ;< NOTE: ENTER  S  OR  D  IN  QUOTES). 
’  D’ 


THESE  PLANS  ARE: 


* : NORMAL : ::::::::::: : TIGHTENED: ::::::::::: REDUCED: :* 


AC  1  = 

2 

AC  1  = 

1 

AC  1  = 

0 

RE  1  = 

5 

RE  1  = 

4 

RE  1- 

4 

Nl  = 

80 

Nl  = 

80 

N1  =* 

32 

AC  2= 

6 

AC  2= 

4 

AC  2= 

3 

RE  2= 

7 

RE  2= 

5 

RE  2=* 

6 

N2= 

80 

N2= 

80 

N2= 

32 

DO  YOU  WANT  SCHEME  EVALUATION.  .  ? 

IF  YES  ENTER . 1 

IF  NO  ENTER . £ 

1 

DO  YOU  WANT  A  TABLE  OR  A  GRAPH  FORMAT 
FOR  GRAPH.  .  .  . ENTER:  1 

F OF  TOfl  E.  ...  ENTER: 
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ftC  1=  2 

RE  1=  5 

80 

AC  1= 

RE  1= 
Kl= 

1 

M 

80 

1 

1 

!  sas 

i  ii 

1 

|  II  II 

! 

0 

4 

32 

K  2=  i 

AC  2= 

4 

AC  2= 

3 

RE  2=  7 

RE  2- 

5 

RE  2= 

0 

K2=  80 

«2= 

80 

K2= 

32 

I 


fWEWCE  SAMPLE  NUMBER 

,320£*02  M0SE«02  i<WC*«2  575E+02  .itHHW  7*&*02  831E+02  ?1£*02  100E+«3  lWE+$3  117E*3 

PERCENT  DEFECTHC*  ♦  ♦♦♦♦♦♦♦♦♦ 

.000* 

.005*  * 

.010*  N 

.015+  * 

020*  » 

.025*  » 

.030*  « 

.035*  * 

OW  « 

,0«i5+  « 

.050+  * 

.055+  « 

Oto*  * 

oa*  * 

.070+  « 

.075*  » 

080+  « 

.085+  « 

0W+  » 

.025+  ■ 

.!«♦  * 


Example  5  QAFI.FOR 


WHAT  DO  YOU  WISH  TO  DO"> 

1 - DERIVE  PLANS  TO  MINIMIZE  ASN 

2- DERIVE  PLANS  TO  MINIMIZE  AFI 

3- EVALUATE  MIL-STD-105D  SAMPLING  SCHEME 

4- EXIT  THIS  PROGRAM 


ENTER  CHOICE 
3 

ENTER  INSPECTION  LEVEL  IN  QUOTATION  MARKS 
E.G.  ;  SPECIAL  s  ’  SI ’  ,  ’  S2'  ,  1  S31  ,  ’  S4’ 

GENERAL  :’l’  3’ 

’  S’ 

ENTER  LOT  SIZE: 

2000 

ENTER  AQL  IN  PERCENT.  REMEMBER,  ONLY  A 

STANDARD  AQL  IS  ALLOWABLE 

1.5 

DO  YOU  WANT  SINGLE  <’ S’)  OR  DOUBLE! 'DM 
SAMPLING  PLANS  ;< NOTE I ENTER  S  OR  D  IN  QUOTES) . 
’D’ 


THESE  PLANS  ARE: 


* : NORMAL :::i:::::::::TI GHTENED :::::::::::: REDUCED : : * 


AC  1=  2 

AC  1  = 

1 

AC  1- 

0 

RE  1=  5 

RE  1<= 

4 

RE  1- 

4 

Nl=  80 

Nl*= 

80 

Nl« 

32 

AC  2*  6 

AC  2= 

4 

AC  2- 

3 

RE  2“  7 

RE  2= 

5 

RE  2- 

6 

N2-  80 

N2- 

80 

N2- 

32 

DO  YOU  WANT  SCHEME  EVALUATION..? 

IF  YES  ENTER . 1 

IF  NO  ENTER . 2 

1 

DO  YOU  WANT  A  TABLE  OR  A  GRAPH  FORMAT  ? 

FOR  GRAPH. . . . ENTER:  1 
FOR  TABLE - ENTER:  2 

2 

SPECIFY  THE  NUMBER  OF  FRACTION  DEFECTIVE  VALUES.. 
12 

ENTER  THE  FRACTION  DEFECTIVE  VALUE <S>, 

<PUT  A  COMMA  BETWEEN  VALUES.) . 

.  015,  .  02,  .  025,  .  03,  .  04,  .  05,  .  06,  .  065,  .  07,  .  0fl,  .  09,  .  1 


SCHEME  OPERATING  CHARACTERISTICS 


P :  :  :  :  : 
0.  015 

: : : :P ( A) : : i 
0. 9943 

i  : : : : ASN : : : : s 
59.  80 

:  s  hOQ : : : : s  s : s : 
0. 0145 

:  AF I 

0.  04 

0.  020 

0. 9467 

95.  47 

0.  0180 

0.  10 

0.  025 

0. 8106 

107. 15 

0. 0192 

0.  23 

0.  030 

0. 5785 

115. 85 

0. 0164 

0.  45 

0.  040 

0. 2833 

1 14. 88 

0. 0107 

0.  73 

0.  050 

0. 1364 

107. 39 

0. 0064 

0.  87 

0.  060 

0. 0621 

99.  40 

0. 0035 

0.  94 

0.  065 

0. 0413 

95.  86 

0. 0026 

0.  96 

0.  070 

0.  0274 

92.  75 

0. 0018 

0.97 

0.  080 

0.  0120 

87.  90 

0. 0009 

0.  99 

0.  090 

0. 0053 

84.  67 

0. 0005 

0.  99 

0.  100 

0.  0023 

82.  65 

0. 0002 

1. 00 

LIMIT 

NUMBER  FOR 

REDUCED  INSPECTION  IS« 

WHAT  DO  YOU  WISH  TO  DO? 

1 — DERIVE  PLANS  TO  MINIMIZE  ASN 

2—  DERI VE  PLANS  TO  MINIMIZE  AFI 

3—  EVALUATE  MIL-STD-105D  SAMPLING  SCHEME 

4— EXIT  THIS  PROGRAM 


ENTER  CHOICE 
4 

FORTRAN  STOP 


APPENDIX  B 
Program  Listin 


Program  QMAIU.FOR 


"  '•  1  -■  i  '  Hmuf**  I-  »«•«*«*»*«««•««  tnt»«f  i'  t  n  «*#«(,«»  »#p 

3CC.  CONTROL!  INC  PROGRAM  (TIP  TFIE  GUAI  I  TV  CMN'ROI 

POO--’  (  INSPECTION  1 '.AMP L  ING  SOI  ’’WART  PACKAOI 

i  )()('<: 

0006  i  OH.  RICHARD  S.  t.  EA  VENWOP  I  H 

000/  c  df.par imem  or  industrial  and  systems  engineering 

0000  L  UNIVERSITY  Of  FLORIDA 

0007  C  GAINFSVII  I  I  ,  FLORIDA  32611 

00  J  1  c 

0017  BYTE  OUff  il  (H) 

0O10  COMMON /BL  K  1  /N2S,  N2I 

OOl  A  <  0NM0N/BI  K2/PS-  PL 

0015  COMMON/BI  K3/N1 

OOI/.  COMMON /Bl  K4 /ALPHA,  BE  LA 

0(0/  COMMON/Bt  K5/P0,  PI 

0010  COMMON/BL K6/C 1 ,  C2 

0017  COMMON/BL  K7/SUMI.0G  (  4000  ) 

C’OPO  COMMON/BL  KO/N 


OOP  1 

COMMON/BL  KV/C2MAX,  Cl  MAXI  15) 

Op,.’  ■ 

CONMON/Bl  K 1 O/NS,  NL 

no  >  . 

COMMON/ BL El  1 /AHN,  ASNMAX 

OO,  -'I 

C0MM0N/U1  K12/0UTFIL 

OOP', 

C 

00P6 

NUM  0 

OOP/ 

1 0 

NUM“NUM+ 1 

OOPO 

WRITE  (5,  15) 

00P7 

1  5 

FORMAT (//////,  23<  '  ),  'WHAT  DO  YOU  WISH  TO  DO?') 

OOOO 

WRITE (5, 20) 

0031 

?0 

FORMAT!//, 15( '  '),  '1-DERIVE  PLANS  TO  MINIMIZE  ASN ' ) 

003? 

WRITE (5,  25) 

0033 

2  5 

FORMAT! /, 15< '  '),  '2-DERIVE  PLANS  TO  MINIMIZE  AFI ' ) 

0034 

WRITE! 5, 30) 

0035 

30 

FORMAT! /. 15! '  '),  '3-EVALUATE  MIL-STD-105D  SAMPLING 

0036 

WRITE! 5,  35) 

0037 

35 

FORMAT! /,  15!  '  '),  '4-EXIT  THIS  PROGRAM') 

0030 

WRITE!  5,  40) 

0037 

40 

FORMAT! /////,  '  ENTER  CHOICE') 

0040 

READ!  5,  45)  ICH 

0041 

45 

FORMAT! I I ) 

004? 

c 

0043 

IF  ((ICH  . NE.  1). AND. (ICH  . NE.  2). AND. (ICH  . NE.  3). 

0044 

% ( ICH  . NE.  4) )  THEN 

0045 

WRITE  (5,  50) 

0046 

50 

FORMAT!/,  15!/),  '  YOU  MUST  ENTER  1,2,3,  OR  4  '  ) 

O  04  7 

GO  TO  10 

0040 

FND1F 

004  7 

c 

0050 

IE!  (ICH  .  NL.  4)  .AND.  (NUM  .EG.  1))  THIN 

0051 

WRITE (5. *)'  WHAT  DO  YOU  WISH  TO  CALL  YOUR  OUTPUT 

0052 

READ (5, 55)  OUTFIL 

0053 

55 

FORMAT ( 10A1 ) 

0054 

CAL  L  ASSIGN(  1,  OUTFIl  ) 

005  5 

f-  ND  I  F 

0056 

c 

0057 

GO  TO  (60-65,  70,75),  ICH 

UNA  1  NT.MA  I  M 

0050 

i 

0057 

> 

CAI  1  tA:.N 

0060 

GO  10  10 

0061 

t'j  \ " 

CAI  1.  OAI  1 

0062 

GO  ro  io 

0063 

/o 

CAI  L  GMI1 

0064 

go  in  j 0 

0065 

/5 

ST(  IP 

0066 

F  N!) 
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Program  Q ASN. FOR 


OU'  ■  i 
Cm  H  • 

000:5 
0004 
ooo:- 
0006 
0007 
0008 
0007 
0010 
001  1 
ooi  7 
ooi  :? 

001  'I 

oo  1  o 
00  1  6 
00  1  / 
00 1  8 
0017 
0020 

002  1 
nfpn 

0028 
0024 
0025 
0026 
0027 
0020 
0027 
0030 
0031 
0032 
0033 
0034 
0035 
0036 
0037 
0038 
0037 
0040 
0041 
0042 
0043 
004  4 
0045 
0046 
004  7 
0043 
0047 
0050 
0051 
0052 
0053 
0054 
0055 
0056 
0057 


•a  "".r-  ■:  »’  C<MI****»*»»IHHt**»(.»»**»(H**##«(HI*#**e 

cvo.v  it-  rum  hoi  ouuout  sampling  program  it:  ana:  yse 

lMlUDU.  '  AMv’l  INC  PLANS.  ABN  (POT  AND  ASNMAX. 

C  RINUMIAI  AND  POISON  PRUBAL I  TV  0  I STR I  DU r I QNS  USED. 

( 

C  PROGRAMED  BY  R.  WAREN  RANGARaJAN 

C  DEPARTMENT  OF  INDUSTRIAL  AND  SYSTEMS  ENGINEERING 

C  UNIVERSITY  OF  FLORIDA 

C  GAINESVILLE,  FLORIDA  32611 

(:»»*««-**##«***#*»«#«»*»##**#*##****nt»#**«**(Ht*«#»#«»*ti 

r 

SUBROUTINE  GA'JN 

DOUBLE  PRECISION  SUMLUG 

INTEGER  C,  Cl,  C2,  Cl  MIN,  C2MIN,  PI  .  R1  1 

BYTE  0UFK1I  (8) 

( 

COMMON/BI  K1/N23,  N2L 
C  OMMON/BI  K2/PS,  PL 
COMMON/BI  K3/N1 
COMMON/Bt  K4 /ALPHA-  BETA 
COMMON/BI  K5/P0.  P  1 
COMMON/ Ul  K6/C  1 ,  C2 
COMMON / B I  K 7 /SUML OG ( 4 000 ) 

COMMON/131. KU/N 

COMMON/BL K7/C2MAX,  C1MAX< 15) 

COMMON/BL. K  10/NS,  NL 
C0MM0N/BLK1 1 /ASN,  ASNMAX 
COMMON/BL K 1 2/OUTF I L 

C  «■**  -a-***  **  ***********************************  *********  * 

C  BEGINNING  INITIALIZATION 

C  ******************  ************************************* 

N--0 

C2=1000 
ASNM I N= 1 5000. 

C=- 1 

C  ******  **************************************  ********** 

C  INPUT  FORMAT 

C *******  **********  ************************************* 

10  WRITE  (5,  15) 

15  FORMAT  <///'  CODES  R3R  SELECTING  APPR.  PROB.  DIST.  '/ 
1  15X,  'BINOMIAL  ',  12X,  '  =  1 
2/15X.  'POISSON',  13X,  '=2') 

READ  ( 5,  *)  K 

IF(K.  GT.  2.  OR.  K.  LT.  1  )  GOTO  10 
20  WRITE (5, 25) 

25  FORMAT ( 1  OX ,  ' SELECT '/ 1 6X ,  'SAMPLE  PLANS  ONLY  =1' 

1/16X,  'ASN  VALUES  ONLY  -2' 

2/16X,  'OR  BOTH  -3 ' ) 

READ  (5,  *)  KOPT 

IF  <  KOPT .  GT . 3. OR. KOPT. LT.  1 )  GOTO  20 
WRITE  ( 5, 30) 

30  FORMAT ( 1  OX  ,  'INPUT  ALPHA  ') 

READ  (5,*)  ALPHA 
WRITE  (5, 35) 

35  FORMAT (10X,  'INPUT  BETA  ') 

READ  (5,  *)  BETA 
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U iV  ,N 


0(  >'  i!  i 

oor.9 
0060 
006  J 
006? 
0060 
()(>6/l 
0060 
0066 
0067 
00613 
0069 
0070 
0071 
0072 
00  713 
007*1 
0070 
0076 
0077 
00713 
0079 
0080 
0081 
0082 
0083 
0084 
0085 
0086 
0087 
0088 
0089 
"0090 
0091 
0092 
0093 
0094 
0095 
0096 
0097 
0098 
0099 
OlOO 
OlOl 
0102 
0103 
0104 
0105 
0106 
0107 
0108 
0109 
01  10 
01  1  1 
0112 
0113 
0114 


( 


WR]  It  (5.40) 

40  FORMAT  (1  OX,  '  JNf'Ut 
ROAD  (5,  *)  lJ() 
WRITE  (5,  45) 

45  FORMAT (1  OX,  'INRUT 
READ  (  5,  *  )  PI 
WRITE! 5, 50) 

50  FORMAT!  5X. 'INPUT 
1 '  IF  NO 


PO  '  t 


PI  '  ) 


A  SEED  FOR  SINGLE  SAMPLING  NO.  '// 

SEED  AVAILABLE  ENTER  ZERO  AS  THE  SEED  VALUE') 


READ  !  5,  *  )  NS 
WRITE! 5, 55) 

55  FORMAT!  5X,  'INPUT  A  VALUE  FOR 
1 '  IF  R 1 -C2  THEN  THE 

2  '  IF  R  1  >C2  THEN  THE 

3  '  IT  Rl! C2  THEN  THE 

READ!  5,  *)  Rl  1 


(R1-C2)  '// 

VALUE  WOULD  BE  0'// 

VALUE  WOULD  BE  A  POSITIVE  NO.  '// 
VALUE  WOULD  BE  A  NEGATIVE  NO.  ') 


WR I TE  ( 1 ■ 60 ) 

60  FORMAT!  '1  '.  ///10X,  'DEPT.  OF  ISE  ' 

1/, 10X,  UNIVERSITY  OF  FLORIDA  '/ 

2/5X,  5!'*'),  'DOUBLE  SAMPLING  SYSTEM  TO  MINIMIZE  ASN ' , 

35!  ' * '  ) ,  2X,  /  ) 

WRITE  !  5,  65  )  ALPHA,  BETA.  PO,  P  1 
WRITE  (1.65)  ALPHA,  BETA,  PO,  PI 
65  FORMAT! //lOX.  'ALPHA  =',F6.  4,  5X,  'BETA  =  ',F6.  4. 

1/10X,  'PO  =  '  i  F6.  4,  8X  ,  'PI  =■  ' ,  F6.  4  > 

WRITE  (5.  70)  R  1 1 
WRITE!  1,70)  R  1 1 

70  FORMAT ( /5X.  'REJECTION  NO.  OF  FIRST  SAMPLE  (Rl)  =  C2+('.I3,  ')') 
MCI- 10. 0/ (PI /PO) 

75  C=C+1 

C  **************  **************************************** 

C  SINGLE  SAMPLING  PROCEDURE  BEGINS 

C ****************************************************** 

80  NS“NS+ 1 

C ********* ********************************************* 

C  COMPUTATION  OF  LOWER  BOUND  OF  SINGLE  SAMPLING  PLAN 

C  ****************************************************** 

IF ( K. EQ. 1 )  CALL  PROBS 1 < NS. P 1 . C , BXLEC, N ) 

IF(K.EQ.  2)  CALL  PR0BS2 < NS,  P 1  -  C,  BXLEC,  N ) 

IF(BXLEC.  GT.  BETA)  GOTO  80 
NLT  =NS- 1 
NL-MAXO ( 1, NLT) 

C ****************************************************** 

C  COMPUTATION  OF  UPPER  BOUND  OF  SINGLE  SAMPLING  PLAN 

C ******************************** *********  ************* 

85  NL=NL+1 

IF ( K. EG. 1 )  CALL  PR0BS1 <NL, PO, C, BXLEC ) 

IF  <  K.  EG.  2 )  CALL  PR0BS2!NL,  PO,  C.,  BXLEC  ) 

IF(BXLEC. GE. ( 1- ALPHA) )  GOTO  85 
NL=NL-1 

C *************************** *************************** 

C  TEST  FOR  FEASIBILITY 

C  ****************************************************** 

IF ( NS. GT. NL )  GOTO  75 
WRITE  (5,90) 
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■itV  ,N 


0113 
O!  J6 
-  011/ 
01  1(J 
01  19 
0120 
0121 
0122 
0123 
0124 
0125 
0126 
0J  27 
0128 
0129 
0130 
0131 
01  32 
0133 
0134 
0133 
O  1  36 
0137 
0138 
0139 
0140 
0141 
0142 
0143 
0144 
0145 
0146 
0147 
0148 
.  0149 
0150 
0151 
0152 
0153 
0154 
0155 
0136 
0137 
0158 
0159 
0160 
0161 
0162 
0163 
0164 
0165 
0166 
0167 
0168 
0169 
0170 
0171 


WE  J  !  I  '  1  /u  ; 

F0RNA1  •  .  /<;x  »• .  •  NGt  C  SAMPl  '  N<  pi  AM  t 
WH  1  rt:  <  5-  93  )  c  .  M!  ■■  NL 
WR 1  7  f  n .  93  >  C,  MS.  NL 

EDEMA T  <  10X,  'ACCEPTANCE  NO.CC)  =  ,12 

1. /10X,  'LOWER  BCJC/ND  ON  N  (NS)  ^',14 

2,  /10X,  'UPPER  BOUND  ON  N  ( NL )  =  ',I4> 

C********************************  ***##*««««#«*#****#*»# 

C  COMPUTATION  OF  DOUBLE  SAMPLING  PLAN  BEGINS  FOR  EACH  VALUE  OF  C2 

(****#*♦**»****««*«*«**♦*#***#*#♦**#***#*♦#**«***«#*#«* 

IFCC.LT.C2)  MC-C+MC1-1 
C2-C 
C 

Rl  CL'OCll 

C 

DO  140  Kl  1  ,  C2 
C  1  -  K  1  1 

C *»«•**•»*******«  **********  #**■**  #■****•»■ -a-****  *•»•**■«■  * »**»##*« 

C  CALL  SUBROUTINE  TO  COMPUTE  THE  FIRST  SAMPLE  NUMBER 

Cim*##****##**********#*****#**#**#**!******#********* 

CALI.  1RY1  (  NTRY ,  C  1.  P  1,  NS,  BFTA,  K> 

Nl-  NTRY 

IFCNTEY. 07.  NS)  GOTO  J 45 
C 
C 

WRITEC  5,  100) 

WRITE ( 1 , lOO) 

1O0  F0RMATC/10X,  'DOUBLE  SAMPLING  PLANS',/) 

WRITEC  5,  105)  Cl,  C2 
WRITEC 1 ,  105)  Cl.  C2 

105  FORMAT  C/10X,  'FOR  C  1  =■  '  »  12,  2X,  'C2~  ' ,  12,  /  /  ) 

C 

NTEMP=N1 

IFCKOPT. EQ. 1 )  WRITEC 1,110) 

IFCKOPT. EQ. 1 )  WRITEC5. 110) 

IFCKOPT. EQ. 2)  WRITEC1.112) 

IFCKOPT. EQ. 2)  WRITEC5,  112) 

IFCKOPT. EQ. 3)  WRITEC5,  115) 

IFCKOPT.  EQ.  3)  WRITEC1.U5) 

110  FORMAT  C11X,  '  (Nl )  ', 3X,  '  C  N2S )  C  N2  <  C N2L )  ' . 

1  6X ,  'PS',  lOX,  'PL'//) 

112  F0RMATC10X,  'NH',3X,  'N2S',4X,  'ASNMAX',6X,  '  ASN  ' ,  /  /  ) 

115  FORMAT  C  1  OX ,  '  CN1  )  ',  3X,  '  CN2S)  ',  3X,  '  C  N2L  )  ',  4X, 

1  'ASNMAX 6X, 'ASN'/) 

NTEMP 1 =NS 

C  ***************************************  *************** 

C  COMPUTATION  OF  SECOND  SAMPLE  FOR  EACH  VALUE  OF  FIRST  SAMPLE 

C#«***»#**********#******#««*»#**#**»*»t»fHHHHt********* 
ASN=FLOATCNS>*10 
DO  135  IZ=NTEMP,  NTEMP 1 
I  =  IZ 

IF  ( C NTEMP 1  — I )  .  LE.  CC2*1.5)>  GOTO  135 
C*»**##***«**«***************************************** 

C  CALL  SUBROUTINE  TO  COMPUTE  SECOND  SAMPLE 

C it**********************************  -****««**«********* 

CALL  TRY2CNS, NL, K,  I, Rl  ) 

IFCKOPT. NE. 1 )  GOTO  125 
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1  2u 


'  ■  ',4 


ro.  6,  4\,  fb.  6' 


(MSN 


0172 
0173 
0174 
017  3 
01  76 
0177 
01  7B 
0179 
01  BO 
0181 
0182 
0183 
0184 
0185 
0186 
0187 
0180 
0189 
0190 
0191 
019? 
0193 
0194 
0195 
0196 
0197 
0198 
0199 
0200 
0201 
0202 
0203 
0204 
0205 
0206 


C 

c 

c 


WRJ  IE<  1  .  120  ■  1  N2S,  N 2I.PS.PL 

WR]1L'(5.  120)  I.  N2S,  N2L  .  PS,  PL 
FORMAT ( 1  OX.  14,  5X,  14  <  N2 

GOT  O  135 


TEST  FOR  FEASIBILITY 

***«*#♦*##»#*« 


125  I F ( N2S. GT - N2L )  GOTO  135 

IT <N25. LT. (C2-C1 ) . OR. I. LE. C2>  GOTO  135 
ASNTEM=ASN 

C  ************  »*#***#****************#*«*#«-********#*#** 

c  CALL  SUBROUTINE  TO  COMPUTE  ASN(PO)  AND  ASNMAX  VALUES 

C*####***##»***»»* «#**«****»**«***««***«# ********■•»****■•& 


CALL  ASNN ( MC ,  NS, K,  I, KQPT,  C1MIN, C2MIN,  N1MIN,  N2MIN,  ASNMIN) 
IT (KOPT. NE. 3)  GOTO  135 
WR IT  E (  1 .  130)1,  N2S,  N2L.  ASNMAX , ASN 
WRITE  (5,  130)1,  N2S,  N2L,  ASNMAX,  ASN 
130  FORMAT  (  10X,  13,  5X,  13,  5X,  1 3,  4X  ,  2  <  F8.  4 ,  3X  )  ) 

IF ( ASN. GT. ASNTEM )  GOTO  140 
135  CONTINUE 
C 

140  CONTINUE 

C 

145  IF (C. LT. MC )  GOTO  75 

I F ( KOPT  .FO.  1)  GO  TO  153 

WRITE ( 1,  150)  ASNMIN,  N1MIN,  N2MIN,  C1MIN,  C2MIN 
WR I TE ( 5,  150)  ASNMIN,  N1MIN,  N2MIN,  C1MIN,  C2MIN 
150  FORMAT(////, 10X, 'GLOBAL  MINIMUM  ASN ( PO ) = ' , F8. 2, / / 

UOX,  'CORRESPONDING  N1  =  ',  15// 

210X,  'CORRESPONDING  N2S  -',15// 

310X,  'CORRESPONDING  Cl  =',  12// 

410X,  'CORRESPONDING  C2  =', 12) 

C 

155  RETURN 
END 


0001 
0002 
0003 
0004 
0005 
0006 
0007 
0008 
0009 
0010 
001  1 
0012 
0013 
0014 
0015 
0016 
0017 
0018 
00  1  7 
0020 
0021 
002? 
0023 
OOP  4 


SUBROUTINE  IRY1 ( NT  R Y ,  Cl,  P,  NL,  BETA,  K) 

C  *»■#■«  -tnunt -turn-****#  *■**■***•*•**•***»  *»***-li-****-fr-****#i«-******** 

C  THIS  SUBROUTINE  COMPUTES  FIRST  SAMPLE  NUMBER  OF  DOUBLE 

C  SAMPLING  PLAN  BY  AN  INTEGER  FORM  OF  BISECTION  METHOD 

C  ***************************  ♦♦■a-*********-*-****#*-**'**#*** 

INTEGER  Cl 
C 

NLARGE=Nl 

NSMALI.-O 

C 

10  NT R Y= ( NSMALl  ENLARGE )/ 2. O 

C  **■**«  *«■*********«•*•»**«■*********«■«•*«■**#*«■«■«•*■«■«■*«■*■»■***#* 

C  CALL  APPROPRIATE  PROB A I L I TY  SUBROUTINE  FOR  PROB.  CALCULATIONS 

C  ■*•»***•***•*********■■***•#■***■***■*****«•****«****■*•*■*■*#*■»■**  *«* 

15  IF(K.EQ.l)  CALL  PROBS 1 ( NT R Y ,  P ,  C 1 ,  B XLEC  ) 

IF ( K. EQ. 2 )  CALL  PR0BS2 ( NTRY, P ,  C 1 ,  BXLFC  ) 

IF1BXLEC. LE. BETA)  GOTO  20 

NSMALL=NTRY 

GOTO  25 

rO  NLARGF-N1RY 

25  IF  ( NSMALL.  Nf-.  (NLARGE  1))  GOTO  10 
N  TRY=NI_ARGF 
RETURN 
f-  Ml) 


Oivj  subroutine  thyzmns.  ml.,  k.  j,  mi  ) 

OOO.  .  >•  »  <*  *  ■  4Mt  tut*  o  #■*«■«•■»•«■  *«■  a  -tut-* -ihmm;  a  ***«*•**■« -tt  tum  *** -a  a  * 

ooc>\;  .  rue;  subroutine  computes  the  second  samp;  e  number  jf 

0004  (HI  DOUBLE  OAMPLINC  NUMBER  Br  AN  I N  CL  GET?  BISECTION 

OOOD  r  METHOD.  SEVERAL  TESTS  ARE  DONE  TO  LOCATE  THE  PARAMETER 

0006  (  A  I  1  rs  TRUE  POSITION. 

OOO/  ( 

0008  INTEGER  C1,C2,R1 

0009  t 

0010  COMMON/ BL K1 /N2S,  N2L 

0011  COMMON/BL  K2/PS,  PL 

0012  COMMON /BLK3/N1 

0013  COMMON/BL K4 /ALPHA.  BETA 

0014  COMMON/BL K5/P0,  PI 

00  ID  COMMON/BL K6/C 1  .  C2 

0016  (, 

OO  17  K  1  -C  1  t- 1 

0018  c*  f  a-*  ■»■;<•**  ihhuhhhuhhhhhhhhhhuhhhhhhhhuhhhhhhj# 

0019  C  SE 1  LIMITS  EOR  COMPUTING  N2S 

0020  C  *'>**«*****•»■■***  **<•«■****♦#«•****#»«■**■»**#«*#******«»#*•»#* 

002  1  NSNALI.  =NS  ~J 

0022  Ml.  ARGE-NSMAL1 

0(123  lUftHfBBBBBfrtHt#*#**#*##*##*##*##**#***#*#*#*****-* 

0024  C  INDEXING  TO  SPECIFY  WHAT  BOUND  ( N2S  OP  N2L ) 

OOZ'D  C  IS  BEING  COMI'UTED 

0(.'26  C*1  IHHUHtaBBtHUHtB-KB*#**#***#**##*##*##**#*#**#****####** 

0027  1  :=} 

0028  C  *#**■»# a*  -a-* *  •#■■*#  *■#■*•**■* **■*»••«••*«•*****•*  *•****•#* 

0029  C  INITIAL.  TEST  AT  EACH  LIMIT 

0030  C  *********  *******•»*•«■*  *******  *************  **  #••*•*•*•*  ■#•#•* 

0031  CALL  PROBD1 (J,  NSMALL,  PI ,  DPROB,  K,  R1 ) 

0032  IF(DPROB. LE.  BETA)  GOTO  40 

0033  C  »tHnntt«»##*#fl#****»#*mm**********#»**»»**t»««»#*** 

0034  C  BISECTION  METHOD 

0035  ClHHHHHUHUHUHUHHUt#*#*#*##***#**#*##*##***************#* 

0036  NLARGE=NL 

»  0037  10  NTR Y= ( NSMALL+NLARGE ) /2.  0 

0038  GOTO  (15. 20), I 

0039  C 

0040  15  CALL  PROBD  1  <  J,  NTRY,  P  1 ,  DPROB ,  K,  R  1  ) 

0041  IF ( DPROB. LE.  BETA)  GOTO  30 

0042  GOTO  25 

0043  20  CALL  PROBD 1 ( J, NTRY, PO,  DPROB , K,  R 1 ) 

0044  IFCDPROB. LT. ( 1-ALPHA) )  GOTO  30 

0045  25  NSMALL=NT  RY 

0046  GOTO  35 

0047  30  NLARGE=NTRY 

0048  35  IF( (NLARGE-NSMALL). GT. 1 )  GOTO  10 

0049  ft****************** -a-***** *■**•#•* 

0050  C  CHECK  THE  INDEX  TO  FIND  WHERE  THE  PROCESS  IS 

0051  C*  «■  *«#-**#****-*-**#■«•*■»■■«•■*■** #*#***■»*  *■#■■*  a-#**** 

0052  GOTO  (40,  45) ,  I 

0053  c*  *******  **■****•**»*******■«■*■«■  *■«■■*  **•**#•  ***•*•«■•«•  *•«•*•*■#•#«■ 

0054  C  CHANGE  THE  INDEX  AFTER  N2S  COMPUTATION 

0055  C  ♦»♦»***♦*«****#  ***** 

0056  40  I  =  I  ■*- 1 

0057 


UW'./iJ 

0057 
0060 
0061 
0062 
0063 
006-1 
0065 
0066 
0067 
00 6B 
006V 
00  70 
007  1 
0072 
0073 
0071 
0075 
0076 


0001  SUBROUTINE  PR0BS2  (  NN,  P>  C,  BXLEC  ) 

0002  C«»»*hhhi-ii  **-*•;<■*-«■**#*■«■  *****•**##•*»•*•«••»  ******  *■#•#■»■  **•* *■»*■»***■« 

0003  C  THIS  SUBROUTINE  COMPUTES  CUMULATIVE  POISON 

OOOI  C  PROBABILITIES 

0005  C***#***-##*#*###*###*#******#*#*#*#**#*******#*###***#-#- 

0006  INTEGER  C 

0007  C 

0008  PP-P*NN 

0007  TERM=1.0 

OOIO  SUM=TERM 

OOl  1  C 

0012  ir(C.EQ.O)  GOTO  15 

0013  DO  10  1=1,  C 

0011  TERM  --TERM*PP/  I 

00 1 5  SUM=SUM+TERM 

0016  10  CONTINUE 

OO  1  V  C 

0018  15  BXLFC=SUM/EXP  U’P ) 

0019  C 

0020  RETURN 

0021  END 


'•  IL'  1  INC  EACH  POSSIBLE:  CASES  TO  LOCATE 

<  I  Hr  LOWER  BOUND  AT  ITS  TRUE  POSITION 

N2S~MAX0  (  0,  N1  ARGE) 

CALL  PR0BD1 (U. NLARGE  PI,  DP ROB ,  K ,  R  1  ) 
PS=DPRGB 
MTEMP=NL.ANGL  1 
NSMALL=MAX0(0,  MTEMP  ) 

NLARGE=N1 
GOTO  10 
15  N2L-NSMA1  1 

CALL  PR0PD1 ( J,  NSMALL,  PO,  DPROB ,  K,  R  1  ) 

pl=-dprob 

CALL  PR0BD1  !  J,  NLARGE,  PO,  DPROB,  K,  R  1  > 

IF ( DPROB. GE. (1 -ALPHA))  M2L=NLARGE 
IF ( DPROB. GE.  < 1  ALPHA) )  PL^DPROB 

C 

50  RETURN 
END 


000 1 

oooy 

0003 
0004 
000 5 
0006 
0007 
0000 
0009 
OOl  O 
001  1 
0012 
OO  1  3 
0014 
OO  1  5 
0016 
0017 
0018 
0019 
0020 
0021 
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SUBROUT  I  Nh.  ASNN(MC,  NS,  K,  N1  1  KLIP  T  ,  C1NIN,  C2MIN,  N1M1N,  N2MIN,  ASNMIN) 

C  THIS  SUUKOU  TINE  COMPUTES  ASN(PO)  VAI  UCS  AND 

C  ASNMAX  VALUES. 

DOUBLE  PRECISION  SUMLOC 
INTEGER  C1MIN, C2MIN 
C 

C0MM0N/BLK1/N2S, N2L 
COMMON/BL  K3/N1 
COMMON/ DL  K4 /ALPHA,  BETA 
COMMON/BL  K5/P0, PI 
COMMON /BL.K6/  1 1 ,  12 
COMMON/BL  K7/SUML0G ( 4000 ) 

COMMON/BI  K8/N 
C0MM0N/BLK1 1 /ASN, ASNMAX 
COMMON/BL  K  1  2/0UTF IL. 

C  •*■#•*•»•*■***•#•#•#  *{*  *#■»****  *  ***■**■#■#■****■***-**#■  *  *  **  **■ * *****  *  **  * 

C  INITIALIZATION 

C  COMPUTE  P*  (MAXIMUM  PROD.  FOR  ASNMAX) 

C*#**-#-**#******  ************************* #*•*********##*•» 

J--  liu 


c 


C 


C 


X  X  X-0. o 

I F ( III G I . O  >  XXX-SUMLGG< 1 1 )  „  .  . 

AK0NST=10.  **(SUMLOG<  12)  <-SUMLOG(Nl  1-12-1  ) -XXX-SUMLOG  (  N1  1-1  1-1  )  ) 
TEMP-1. O/PLOAT ( 12-11) 


AKONST=AKONST«-«-TEMP 
PSTAR=AKONST / ( 1. O+AKONST ) 

IF ( K. EQ. 1 )  CALL  PR0BS1 <N1 1 , PSTAR, 12, BXLEC ) 
IF ( K. EG. 2 )  CALL  PR0BS2 (N1 1 , PSTAR, 12, BXLEC ) 
TEMP=BXLEC 

IF ( K. EQ.  1 )  CALL  PR0BS1 < N 1 1 , PSTAR,  1 1 , BXLEC ) 
IF ( K. EQ. 2 )  CALL  PR0BS2 ( N1 1 , PSTAR,  1 1 . BXLEC ) 
TEMP 1=TEMP-BXLEC 
ASNMAX=FLOAT (Nil) +N2S*TEMP 1 


IF ( K. EQ. 1 )  CALL  PR0BS1 ( N1 1 , PO, 12, BXLEC ) 

IF ( K. EQ. 2 )  CALL  PR0BS2 ( N 1 1 , PO. I 2, BXLEC ) 
TEMP=B  XL  EC 

IF (K. EQ.  1 )  CALL  PR0BS1 (N1 1 , PO,  1 1 ,  BXLEC ) 

IF ( K. EQ. 2 )  CALL  PR0BS2 < N1 1 , PO,  1 1 ,  BXLEC ) 
TEMP2=TEMP -BXLEC 
ASN^FLOAT (Nil) +N2S*TEMP2 

IF( ASNMAX. GT. NS. OR. ASN. GT. ASNMIN)  GOTO  10 

ASNMIN=ASN 

N1MIN=N1 1 

N2M I N=N2S 

C 1 M I N= I  1 

C2MIN=I2 


10  IF ( KOPT. NE. 2 )  GOTO  20 

WRITE! 1,15)  N1 1. N2S, ASNMAX, ASN 
WRITE(5, 15)  N1 1 , N2S, ASNMAX, ASN 
15  FORMAT (  10X.  2(  13,  3X>.  2<F8.  4.  3X>  ) 

20  RETURN 
FND 


» 


I 


I 


I 


| 


» 
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Program  QAFI.FOR 
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C. 

C  ***»*#**•»*****#•****••*  ***********************  •*■******#*«■* 

C  QUALITY  CONTROL  PROGRAM  TO  DERIVE  DOUDLE  SAMPLING 

C  PLANS  TO  MINIMIZE  AVERAGE  FRACTION  INSPECTED. 

C 

C  PROGRAMMED  BY  JO  ELLEN  WALKER 

C  DEPARTMENT  OF  INDUSTRIAL  AND  SYSTEMS  ENGINEERING 

C  UNIVERSITY  OF  FLORIDA 

C  GAINESVILLE,  FLORIDA  32611 

C  **•»•■***■»**#**####**#****#  ###**#*#**#*******•******#**«•«•* 

c 

SUBROUTINE  OAF I 

INTEGER  Cl, C2, C. CSTAR,  C2M1 

INTEGER  Rl, R1M1 

DOUBLE  PRECISION  SUMLOG 

BYTE  OUTF IL ( 8 ) 

COMMON/BLK4/ ALPHA,  BETA 
C0MM0N/BLK6/C 1 ,  C2 
C0MM0N/BLK7/SUML0G ( 4000 ) 

C0MM0N/BLK8/N 

C0MM0N/BLK12/0UTFIL 

C  *************«**************«  *************************** 

C  INPUT  PARAMETERS 

C ****w************************************************* 
KINDEX=1 
WR I TE (  5,  10) 

10  FORMAT ( 10X, 'ENTER  VALUE  OF  ALPHA') 

READ (5, *> ALPHA 
WRITE  <5,  15) 

15  FORMAT ( 10X,  'ENTER  VALUE  OF  P0'> 

READ  (5,  *>P0 
WRITE! 5, 20) 

20  FORMAT! 1 OX,  'ENTER  AQQL  VALUE') 

READ!  5,  * )  AOGL 
WRITE !  5,  25) 

25  FORMAT! 1 OX. 'ENTER  LOT  SIZE') 

READ! 5,  #)NNN 
C 
C 
C 

WRITE! 1. 28) 

28  FORMAT! '1 ', ///10X,  'DEPT.  OF  ISE  ' 

%/,  10X,  'UNIVERSITY  OF  FLORIDA  '/ 

*/5X, 5!'*').  'DOUBLE  SAMPLING  SYSTEM  TO  MINIMIZE  AFI  ' 
*5!  '*' ).  2X,  /) 

WRITE! 5. 30) ALPHA 
WRITE! 1,  30 > ALPHA 

30  FORMAT! 10X,  ' ALPHA= ' . 2X, F8. 6 ) 

WRITE !  5,  35)P0 
WRITE! 1, 35)P0 

35  FORMAT!  10X,  'P0=  ',  2X,  FS.  6 ) 

WR I TE ( 5, 40 ) AOGL 
WRITE! 1, 40) AOGL 

40  FORMAT!  10X,  ' AOQL= ',  2X,  F8.  6 ) 

WRITE ! 5, 45 )NNN 
WRITE! 1, 45 )NNN 

45  FORMAT!  10X,  'N='.2X,  16) 
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U/V  ] 


0058 
0059 
0060 
0061 
0062 
0063 
0064 
0065 
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0067 
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0079 
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0083 
0084 
0085 
0086 
0087 
0088 
0089 
0090 
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0093 
0094 
0095 
0096 
0097 
0098 
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0108 
0109 
0110 
0111 
0112 
01 13 
0114 


C  #*■*****#*#***«  **************************************** 

C  COMPUTE  SINGLE  SAMPLING  PLAN 

C  ****************************************************** 

c 

C  INITIALIZATION 

C 

C  ****************************************************** 

NS=  1 

C=0 

N=0 

AF  1  =  1.  DO 

C 

WRITE (5.  50) 

WRITE ( 1 . 50) 

50  FORMAT ( /// 10X,  '***  SINGLE  SAMPLING  PLAN  *** ' ) 
C****************************************************** 

C  FIND  NS. C  COMBO  THAT  SATISFIES  L(PO)  G.  T.  1 -ALPHA 

C ****************************** ************************ 

55  CALL  PR0BS1 (NS.  PO,  C.  BXLEC ) 

IF ( BXLEC. LT. (1. DO-ALPHA) >C=C+1 
IF  (BXLEC.  LT.  <  1.  DO-ALPHA)  )NS=*C+1 
C****************************************************** 

C  SEARCH  TO  FIND  MIN  NS  VALUE  SUCH  THAT  AOQL  L.  T.  AOQL* 

C****************************************************** 

60  CALL  SEARCH < NNN.  C,  NS,  SAOQ) 

IF ( SAOQ. LE. AOQL) GOTO  75 
NSTEMP=NNN 

CALL  SEARCH < NNN,  C.  NSTEKP.  SAOQ) 

I F ( SAOQ.  CT. AOQL ) C=C  +  1 
C 
C 

I F ( SAOQ . CT . AOQL ) NS«C * 1 
C 
C 

BL=NS 

BH»NNN 

65  NSTEMP«I ID I NT ( ( BL+BH ) /2. DO ) 

CALL  SEARCH ( NNN,  C ,  NSTEMP , SAOQ ) 

IF ( SAOQ. LE. AOQL ) BH-NSTEMP 
C 

IF  < SAOQ. LE. AOQL ) AOQ-SAOQ 
C 

IF < SAOQ.  GT.  AOQL ) BL “NSTEMP 
IF< (BH-BL).EQ. 1. DO) GOTO  70 
GOTO  65 
70  NS=NSTEMP 

IF ( SAOQ.  GT. AOQL ) NS=BH 


C  CHECK  THAT  NS, C  COMBO  STILL  SATISFIES  L(PO)  CONSTRAINT 

C****************************************************** 

75  CALL  PR0BS1  (NS,  PO,  C,  BXLEC) 

IF ( BXLEC. GE. ( 1. DO-ALPHA) ) GOTO  80 
C=C+1 
C 

NS=C+1 

C 

GOTO  60 


4? 


0115  C  ************************  *****************  ************* 

0116  C  COMPUTE  API 

0117  C****************************»************************* 

0118  80  ATIPO=NS*BXLEC+NNN*< 1.  DO-BXLEt) 

0119  AFIPO=ATIPO/NNN 

0120  C 

0121  NSTAR=NS 

0122  C 

0123  WRITE ( 5*  85  )NS 

0124  WRITEd.  85) NS 

0125  85  FORMAT ( //10X<  'NS='. 13) 

0126  WRITE ( 5<  90 ) C 

0127  WR I TE ( 1 .  90 ) C 

0128  90  FORMAT ( 10X<  'C=',I2) 

0129  WRITE( 5.  95) AFIPO 

0130  WR ITE( 1 .  95) AF IPO 

0131  95  FORMATt 10X.  'AFI (P0)= ' ,  F8. 6) 

0132  C 

0133  100  CSTAR=C 

0134  C ********************************************* ********* 

0135  C  DOUBLE  SAMPLING 

0136  C ****************************** ************************ 

0137  C 

0138  WRITE ( 5. 105) 

0139  WRITE(1» 105) 

0140  105  FORMAT ( / / / 1 OX .  '**  DOUBLE  SAMPLING  PLANS  **  '  ' 

0141  PLAN=1 . DO 

0142  C2-CSTAR 

0143  110  R1=C2+1 

0144  C2M1=C2— 1 

0145  C1P0=C1+1 

0146  R1M1— Rl— 1 

0147  DDATI=*=NNN 

0148  TTCMIN=1. 

0149  C 

0150  C1=0 

0151  JJ=0 

.0152  C 

0153  WRITEI5.  115) 

0154  WR I TE  < 1 .  115) 

0155  115  FORMAT  <  ///  10X*  'C1',6X,  'C2'.7X,  'Nl'.BX.  'N2 ' .  9X»  '  AF  I  ' .  /  /  ) 

0156  C****************************************************** 

0157  C  CALCULATE  FIRST  SAMPLE  NUMBER 

0158  C 

0159  C  FROM  RESULTS  OF  PREVIOUS  RUNS.  IT  WAS  FOUND  THAT  N1  IS  NOT  LESS 

0160  C  THAN  NSTAR/B.  THUS.  THE  INITIAL  VALUE  OF  N1  IS  SET  ACCORDINGLY. 

0161  C****************************************************** 

0162  C 

0163  120  DO  165  LL=INT(NSTAR/8) .  NSTAR 

0164  N1=LL 

0165  IF <N1 .  LT.  C2)N1=»R1M1 

0166  C ****************************************** ************ 

0167  C  CHECK  BtNl.PO.  C2)  G.  T.  1-ALPHA  CONSTRAINT 

0168  C****************************************************** 

0169  125  CALL  PR0BS1 (N1 ,  PO,  C2.  BXLEC ) 

0170  IF < BXLEC. LT.  d. -ALPHA) ) GOTO  175 

0171  C****************************************************** 


C  CALCULATE  SECOND  SAMPLE 

C  ***************  *********  ***********■«■#*■»****•»##*#  **#«•** 

N2=NSTAR-N1 

C**********##*####*#*******#**#**#*******#************* 

C  CHECK  THAT  DOUBLE  PROBABILITY  Q.J.  i-ALPHA 

C-*##*##**####*-*###***##*****##**##*****#-#***#-**##**  #«  *  H 
130  CALL  PR0BD1 (Nl. N2, PO, DPROB, KINDEX, R1 ) 

IF (DPROB. GE. ( 1. -ALPHA) JGQTQ  135 
IF ( C 1 . EQ. C2M1 >  GOTO  175 
C1=C1+1 
GOTO  130 

C****##**#*****#********#***##****##***********#**#**** 

C  FIND  N2  VALUE  SATISFYING  AOQL  L.  T.  AOQL* 

C ****************************************************** 

135  CALL  SRCH2(NNN>  Nl» N2» A0Q1 ) 

IF(A0Q1.LE. AOQL) GOTO  150 

C****************************************************** 

C  N2  WILL  NOT  BE  LESS  THAN  Nl*9,  THE  INITIAL  LOWER  BOUND  ON  N2 

C ****************************************************** 

N2TEMP=N1*9 

C 

CALL  SRCH2<NNN,  Nl.  N2TEMP,  A0Q1 ) 

IF < A0Q1. GT.  AOQL) GOTO  165 
C 

BL=N2 

BH=N1*9 

140  N2TEMP=INT< <BL+BH>/2> 

CALL  SRCH2(NNN.  Nl,  N2TEMP.  AQQ1 ) 

IF ( A0Q1 . LE- AOQL ) BH=N2TEMP 
I F ( AOQ 1 . LE. AOQL ) FAOQ=AOQ 1 
C 

IF ( A0Q1 . GT. AOQL ) BL=N2TEMP 
IF  < ( BH-BL ) . EQ. 1 .  )GOTO  145 
GOTO  140 
145  N2=N2TEMP 

IF ( A0Q1 . GT. AOQL ) N2=BH 

C****************************************************** 

C  CHECK  THAT  BINOMIAL  PROBABILITIES  ARE  G.  T.  1-ALPHA 

C*******#***#*#*###**#*»##»#**»##»#********«*****«****« 

CALL  PROBDKNl,  N2,  PO,  DPROB,  KINDEX.  R1  ) 

IF(DPROB.LT. ( 1. -ALPHA) ) GOTO  165 
150  CALL  PR0BSKN1.  PO,  Cl.  BXLEC) 

PA2=tDPROB— BXLEC 

C ****************************************************** 

C  COMPUTE  ATI 

C****************************************************** 
DATI*N1*DPR0B+N2*PA2+NNN#< 1.  -DPROB ) 
C****************************************************** 

C  IF  THE  ATI  INCREASES,  CONTINUE  FOR  5  ADDITIONAL  INCREASING 

C  ITERATIONS.  THEN  INCREMENT  Cl  AND  CONTINUE. 

C ****************************************************** 

IF< JJ. EQ. 4) GOTO  155 

IF (DAT I. GE. DDATI ) JJ=UJ+1 

I F ( DAT I . GE. DDAT I ) GOTO  165 

DDATI=DATI 

DDAOQ=FAOQ 


UA)  J 


0229 

0230 

0231 

0232 

0233 

0234 

0235 

0236 

0237 

0238 

0239 

0240 

0241 

0242 

0243 

0244 

0245 

0246 

0247 

0248 

0249 

0250 

0251 

0252 

0253 

0254 

0255 

0256 

0257 

0258 

0259 

0260 

0261 

0262 

0263 

0264 

0265 

.0266 

0267 

0268 

0269 

0270 

0271 

0272 

0273 

0274 

0275 

0276 

0277 

0278 

0279 

0280 

0281 

0282 

0283 

0284 

0285 

0286 

0287 

0288 

0289 

0290 

0291 

0292 

0293 

0294 

0295 

0296 

0297 

0298 

•0299 

0300 

0301 

0302 

0303 

0304 

0305 

0306 

0307 


K1=C  1 
K2=C2 
K3=N  1 
K4=N2 
C 

GOTO  165 

C **************************** ***** ************  ********* 
C  MINIMUM  OF  CELL  ( TCMIN )  FOUND 

C****************************************************** 
155  TCMIN=DDATI/NNN 
DDAT I =NNN 

C 

WRITE!  5»  160)K1>  K2<  K3,  K4,  TCMIN 
WRITE! 1,  160 )K1, K2,  K3,  K4,  TCMIN 
160  FORMAT! 10X,  12,  6X,  12,  5X,  14,  6X,  14,  7X.  F8.  6) 

C ************************************************ ****** 
C  IF  MINIMUM  OF  COLUMN  IS  FOUND,  INCREASE  C2 

C ********* ****** ********* ****** ********* ********* ****** 
IF ! TCMIN. GE. TTCMIN  >  GOTO  170 
KK1=K1 
KK2=K2 
KK3=K3 
KK4=K4 

TTCMIN=TCMIN 
TT  AOQ=DDAOQ 


IF  !C1.EQ. C2M1)  GOTO  175 
C1=C1+1 
C 

JJ=0 

N1=KK3 

C 

GOTO  125 
C 

165  CONTINUE 

C****************************************************** 

C  MINIMUM  OF  COLUMN  ! TMIN)  FOUND 

C****************************************************** 

170  TMIN=TTCMIN 
TAOG*=TTAOG 

C ****************************************************** 

C  IF  MINIMUM  SAMPLING  PLAN  FOUND.  STOP 

C****************************************************** 

I F ( TM I N. GE.  PLAN ) GOTO  1 90 
C 

PLAN=TMIN 
KKK1=KK1 
KKK2=KK2 
KKK3=KK3 
KKK4=KK4 
175  C2=C2+1 

C ****************************************************** 

C  NEW  BOUNDS  ON  SAMPLING  PLAN  CALCULATED  FOR  NEW  VALUE  OF 

C****************************************************** 

180  CALL  SEARCH ! NNN,  C2,  NSTAR,  AOQ) 

IF ! AOQ. LE. AOQL ) GOTO  185 
IF ( NSTAR. GT.  NNN) GOTO  175 

NST  AR=NST  AR  + 1 
GOTO  180 


185  CALL 
GOTO 


PROBSi !NSTAR, PO, C2, BXLEC ) 
1 10 


190  WR I TE !  5,  195) 

WRITE! 1,  195) 

195  FORMAT! //10X,  'SAMPLING  PLAN  MINIMUMS') 
WRITE! 5, 200 ) KKK 1 ,  KKK2 
WRITE! 1, 200 ) KKK 1 ,  KKK2 
200  FORMAT! /10X,  'Cl=',  12,  2X,  'C2=',  12) 

WR  I  TE !  5,  205 )  KKK3,  KKK4 
WR I TE ! 1 , 205 ) KKK3,  KKK4 
205  FORMAT!  10X,  'Nl=  ',  13,  2X,  'N2=  13) 

WRITE! 5, 210 ) PLAN 
WRITE! 1, 210JPLAN 

210  FORMAT  (  10X,  'MINIMUM  AFI*',FB.6> 

215  RETURN 
END 
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SUBROUTINE  SEARCH ( NNN,  C,  NS.  AOQ ) 

C  ************** -ft************  *********************  «-«***« 
C  SEARCH  TO  FIND  VALUE  OF  PSTAR  USING  GOLDEN 

C  SECTION  METHOD.  INITIAL  LIMITS  OF  0  AND  1 

C  **■•#*##**###***#•**####•*•#•*#*  *•**##•*■*  *#****#*#  ***■*•#■*  *#**■*# 
INTEGER  C 

DOUBLE  PRECISION  SUMLOG 
C  OMMON / B  LK7 / SUMLOG < 4000 ) 

C0MM0N/BLK8/N 

C 


C 

C 


C 

C 


C 


C 


C 


C 


A 1 =0. DO 
A2=l. DO 
T= 1 . D-3 

R  =  5.  D-1*(DSQRT ( 5.  D0)-1. DO) 
H=A2-A1 


PLFT=A1+(R*R)*H 

PRT=A1+(R*H) 


CALL  PR0BS1 (NS, PLFT, C. BXLEC ) 
ATI=(NS*BXLEC  > +NNN* ( 1. DO-BXLEC ) 
AFI 1=ATI/NNN 
A0Q1=PLFT*( 1.  DO-AFI 1 ) 

CALL  PR0BS1 (NS. PRT,  C, BXLEC) 

AT I  =  ( NS*B  XLEC ) +NNN* ( 1. DO-BXLEC) 
AFI2=ATI /NNN 
A0Q2=PR  T*  ( 1 . DO- AF 1 2 ) 

GOTO  110 


100  CALL  PR0BS1 (NS.  PLFT,  C,  BXLEC) 

ATI=>(NS*BXLEC  )+NNN*(  1.  DO-BXLEC  ) 

AFI1=ATI/NNN 

A0Q1=PLFT* (1.  DO-AFI 1 ) 

GO  TO  110 


105  CALL  PR0BS1 (NS.  PRT.  C. BXLEC) 

AT I - ( NS*B XLEC ) +NNN* ( 1 . DO-B  XLEC ) 
AFI2=ATI /NNN 
A0Q2=PRT* ( 1 . DO-AF I 2 ) 


110  IF ( A0Q1. LT.  A0Q2)  GOTO  115 
A2=PRT 
H— PRT— A 1 

IF(ABS(PRT-PLFT).LE.  T)GOTO  120 


PRT=PLFT 
PLFT=A1+(R*R)*H 
A0Q2=A0Q1 
GO  TO  100 
115  A1=PLFT 

H=A2-PLFT 


IF ( ABS(PRT-PLFT ) .  LE.  T ) GOTO  120 

PLFT=PRT 

PRT=A1+R*H 

A0Q1=A0Q2 
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SUBROUTINE  SRCH2 ( NNN, N1 ,  N2 .  AOG ) 

C 

INTEGER  Cl, C2, R1 
COMMON/BL  K6/C 1 ,  C2 

KINDEX=1 
A 1  =0. 

A2=  1 . 

T=. 0001 

R=. 5* ( DSGRT ( 5.  D0)-1.  ) 

H=A2-A 1 

PLFT=A1-MR*R)*H 

PRT=A1-MR*H> 

C 

C 

CALL  PR0BS1 (Nl,  PLFT,  Cl,  PA1 ) 

CALL  PR0BD1 (Nl,  N2,  PLFT,  DPROB,  KINDEX.  R1 ) 
PA2=DPR0B-PA1 

ATI=DPR0B*N1+PA2*N2+NNN*( 1.  -DPROB  > 

AFI 1=ATI/NNN 
A0Q1=PLFT* ( 1. -AFI 1 > 

CALL  PR0BS1 <Ni ,  PRT,  Cl,  PA1 ) 

CALL  PR0BD1 (Nl, N2,  PRT,  DPROB.  KINDEX, R1 ) 
PA2=DPR0B-PA1 

ATI=DPR0B*N1+PA2*N2+NNN*( 1-  -DPROB) 
AFI2=ATI /NNN 
A0Q2=PRT*( 1.  -AFI 2) 

GOTO  110 
C 
C 

100  CALL  PR0BS1 (Nl, PLFT, Cl, PA1 ) 

CALL  PR0BD1 (Nl,  N2,  PLFT,  DPROB, KINDEX.  R1 ) 
PA2=DPR0B-PA1 

ATI*DPR0B*N1 +PA2*N2+NNN* ( 1 .  -DPROB ) 

AFI 1=ATI/NNN 
A0Q1=PLFT*( 1. —AFI 1 ) 

GOTO  110 

105  CALL  PR0BS1 (Nl,  PRT,  Cl,  PA1 ) 

CALL  PR0BD1  (N1.N2.  PRT,  DPROB,  KINDEX,  R1  > 
PA2=DPR0B-PA1 

ATI®DPR0B*N1+PA2*N2+NNN*( 1.  -DPROB) 

AFI2-ATI/NNN 

A0Q2=PRT* ( 1 .  — AF 12) 

C 

C 

110  IF( A0Q1. LT. A0Q2)C0T0  115 
C 
C 

A2=PRT 

h=prt-ai 

IF ( ABS ( PRT-PLFT ) .  LE.  T ) GOTO  120 
PRT=PLFT 
PLFT=A1+(R*R)*H 
A0Q2— A0Q1 
GOTO  100 
115  A1=PLFT 
C 

C 

H=A2-PLFT 

IF ( ABS ( PRT-PLFT ) .  LE.  T )GOTO  120 
PLFT=PRT 
PRT=A1+R*H 
A0Q1=A0Q2 
GOTO  105 
C 
C 

120  PS= ( PLFT  +PRT  > /2. 

CALL  PR0BS1  (Nl.  PS,  Cl.  PA1  ) 

CALL  PR0BD1 (Nl,  N2,  PS.  DPROB, KINDEX, R1 ) 
PA2=DPR0B— PA1 

AFI® ( DPROB *N1+PA2*N2+NNN*( 1. -DPROB) ) /NNN 
A0Q=PS*(1. -AFI) 

RETURN 

END 
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Program  QMIL.FOR 


O'  >'  , 

ooo ; 

1  > 

ooo; ; 

f» 

0 CO 

ooo: 

(.  • 

0006 

(  B  . 

O  00/ 

C  "< 

0000 

C  B  i> 

0009 

C  **> 

ooio 

C*fr 

OO 1  1 

c  *  «• 

OO 1  2 
0013 
0014 

OO  1  T, 

OO  1  A 

00  l  / 

0!)  !  n 
0019 
002!, 
0021 
0022 
0023 
0024 

c 

O  025 

c 

0026 

c 

0027 

c 

0028 

c 

0029 

c 

0030 

c 

0031 

c 

0032 

0033 

c 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

004  6 

004  7 

0048 

004  9 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 


-  *■  *>  o  "  t:  fl  «  &  tt  1  •«»««*-»»«««  «»«»««««»«««*■«*•*««« 


n  j  i  s  i  u.  i  Of*  *  b 

-  -  - 

'  l  11  “S-  !  M  rKBA!  r  l  VI  PROGRAM  1  Cl  LVALUAl!  ** 

mu  mu  .in  loon  sampi  jnc;  scheme.  <m 

v.  t- 

PKOGRAMMEP  UY  AZMA  I  H.SIDDJQ1  *< 

INDUSTRIAL  AND  SYSTEMS  ENGINEERING  DEPT  ** 

UNIVERSITY  OF  FLOP  I  DA,  GA  I  NESV  I  L  l.  E,  ** 

FLORIDA  32611  »« 


+  *■»■**■*■  +  * 

SU1J ROUTINE  QMIL 
DYTE  OUT  f-  I L  (  8 ) 

CHARACTER  I  *2,  CDL  ,  SP 

INTEGER  FS<  FA.  SA,  AR,  SR,  AC.  AT,  AE,  SC,  ST,  Z.  A,  B 
INI  EGER  NN,  NFS,  NFA,  NAE,  NSA,  NAR,  NAC,  NAT,  TN,  TFS,  TEA,  I  AE,  T 
INTEGER  IAR,  TAG,  TAT,  RN,  RFA , RAE,  RAR.  RAC,  RAT,  FRR,  FT  R, FNR. 

*  ,  STR,  SRR.  RF5,  SI,  S2.  S3,  S4,  S 
COMMON /3CL/L  ,  SP,  S,  I ,  NO,  T,  R,  CDL,  FS.  D,  J,  K.  AOL.  M,  FA,  SA, 

t  AR.  SR,  AC,  AT,  AE  SC,  ST,  2,  NTN ■  NFS,  NFA,  NAE,  NSA,  NAR,  NAC,  NAT  , 

*  TN,  TFS,  TT A,  T  AL ,  TEA.  TAR,  TAC  TAT.  RN,  RFS,  RFA,  RAE,  RAR.  RAC,  RAT 
T.  ,  FRR,  FTR ,  FNR,  SNR,  STR,  SRR,  NNNN 

COMMON /Bl  K1  2/OUT  r TL 


THIS  SECTION  ASKS  THE  USER  TO  ENTER  VALUES  FOR 
THE  NECESSARY  PARAMETERS. 


PRINT  *.  'ENTER  INSPECTION  LEVEL  IN  QUOTATION  MARKS  ' 
PRINT  *,  'E.G.i  SPECIAL  "SI".  "S2",  "S3",  "S4"  ' 
PRINT  *,  '  GENERAL  :  "1",  "2".  "3"  ' 

READ  *,  I 

PRINT  *,  'ENTER  LOT  SIZE:  ' 

READ  *,  L 

PRINT  *.  'ENTER  AQL  IN  PERCENT.  REMEMBER,  ONLY  A' 

PRINT  *,  'STANDARD  AQL  IS  ALLOWABLE' 

READ  *. AQL 

PRINT  *.  'DO  YOU  WANT  SINGLE!  "S' ')  OR  DOUBLE!  "D")  ' 

PRINT  *,  'SAMPLING  PLANS  ;  ! NOTE :  ENTER  S  OR  D  IN  QUOTES). 
READ  *,  SP 
WRITE! 1. 2) 

2  FORMAT!  '  1  '.  ///10X.  'DEPT.  OF  ISE  ', 

1  /.  lOX,  'UNIVERSITY  OF  FLORIDA  ' / /.  5X .  5 !  ' # ' )  , 

2  'SAMPLING  SYSTEM  TO  EVALUATE  M I L-STD - 1 05D ' , 

3  5!  '*•'),  2X.  /  ) 

WR IT E  <  1 . 3 )  I 

3  FORMAT!//,  10X,  'FOR  INSPECTION  LEVEL  ' ,  A3 ) 

WRITE! 1,4)  L 

4  FORMAT!  16X.  'LOT  SIZE"  ',18) 

WRITE!  1,6)  AQL 

6  FORMAT! 16X,  'AQL=  '.F10.5) 

WRITE!  1,8)  SP 

8  FORMAT! 16X,  'SAMPLING  PLAN', A3) 
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Ift  tT. 


I  HI-  :dc  Ft.  t  •> is 


oo?o 
0009 
0060 
0061 
006? 
006? 
0064 
006'.- 
0066 
006/ 
0060 
006'? 
0070 
007  1 
007? 
00711 
0074 
00/0 
0076 
0077 
0070 
0079 
0080 
0081 
008? 
0083 
0084 
0080 
0086 
0087 
0088 
0089 
0090 
0091 
0092 
0093 
0094 
0095 
0096 
0097 
0098 
0099 
OlOO 
0101 
0102 
0103 
0104 
0105 
0106 
0107 
0103 
0109 
0110 
01  1  1 
01  12 
01  13 
Ol  14 


Wi  I  >  I  it 

u!  O  (5,7) 
li  Hi  AT  (  /  /  l  (  X 

v.  i  com 

61.  L  S3 
r  il-L  1  Ml  >1  >, 

Al  L  VA!  Ol  7- 
3  O 

IF (RFA.  I  I  .  I OO ) THEN 
F KR~RFA  i- 1 
END  IF 

IF (RAC.  I  I  -  lOO ) THEN 
RFA-RAC 
FRR-RFAK1 
[NO  IF 

IF (RAT. L  T.  100) THEN 
RF  A--RA'! 

I-  RR=  RFA  i  ? 

END  I  F 

IF (RAF. 1 T. 100) THEN 
RFA=RAI 
I  HR=RF A i 5 
END  I F 

I F ( R AR .  I  T.  100) THEN 
RFA -RAR 
FRR-RFA  *  4 
END  IF 

IF(TFA. LI . 1 00 ) THEN 

FTR=TFA+-1 

ENDIF 

I F  <  T  A  T .  LT.  1  00  )  THEN 

TFA=TAT 

FTR-TFA+-2 

ENDIF 

IF ( T AC . LT. lOOJTHEN 
TFA=T  AC 
FTR=TFA+3 
ENDIF 

IF (TAR. LT. 1 00 ) THEN 

TFA=TAR 

FTR— TFA+4 

ENDIF 

IF  <  TAE. L  T.  1 OO ) THEN 

TFA=TAE 

FTR=TFA+5 

ENDIF 

IF (NFA. LT. 100)THEN 

FNR=NFA+ 1 

ENDIF 

IF (NAT. I T. 100) THEN 

NFA=NA1 

FNR=NFA+-2 

ENDIF 

I F  <  NAC . LT.  1 00 ) THEN 

NFA=NAC 

FNR=NFA  *-3 

CNDIF 


w  i  i 

o  1  1  6 

o  I  J  / 
> )  i  i  i  i 
0  1  1  V 

oi  ;?o 
o  1 2 1 
012^ 
0J  23 
01  24 
0125 
0126 
0127 
012(3 
0127 
0130 
O  1  3  1 
0  I  32 
O  1  33 
0134 
0135 
0136 
0137 
0  1  33 
0139 
0140 
0141 
0142 
0143 
0144 
0145 
0146 
0147 
0148 
0149 
0150 
0151 
0152 
0153 
0154 
0155 
0156 
01  57 
0158 
0159 
0160 
0161 
0162 
0163 
0164 
0165 
0166 
0167 
0168 
0169 
01  70 
0171 


C 

C 

C 

C 

C 

C 

C 

C 

C 


1  I  l  l-IAk  .  I  r  !  ■  •  I  HI  ' 


Nl  A  NAN 
1  NR  NAH  '  /! 
t-.NI)  1 1 

1  1  (  NAt  .11  :  >'■>  :  131  f'J 

NFA-NAI 
F NR -NFA  *  6 
END  I  F 

1 F  <  NSA. L  I .  1 OO) THEN 
NTN=NSA 


SNR-NTN ♦ 1 
ENDIF 

IF ( TSA. L I . 100) THEN 

TN-TSA 

5TR - TN  *  1 

ENDIF 

IF ( SC.  I  I .  1 OO  >  I  HFN 
Z  SC 
3RR  2  t-3 
ENDIF 

IF  (ST.  1  I  .  1  00)  HHF.  N 

7-  ST 

SRR=7*2 

ENDIF 

IF  <  SR.  1.1  .  lOO)  THEN 
Z-“SR 
5RR=7+4 
ENDIF 

I F ( RFS. GF.  L.  OR.  FS.  GE.  L.  OR. NQ. GE. L. DR. RN. GE. L)THEN 


THIS  SECTION  CONTAINS  THE  FORMAT  STATEMENTS 
NEEDED  TO  SET  UP  TABLES  FOR  THE  PLAN  STATISTICS 
AT  THE  DIFFRENT  INSPECTION  LEVELS. 

THE  USER  IS  ASKED  IF  SCHEME  EVALUATION  IS 
DESIRED. IF  SO  THEN  SUBROUTINE  DC  IS  INVOKED. 


WRITE (1. 12) 

WR  I TE  (  5.  12) 

12  FORMAT! 10X,  'USE  100  PERCENT  INSPECT ION ' ,  / ,  'AS  SAMPLE', 
%  '  SIZE  EQUALS  OR  EXCEEDS  LOT  SIZE') 

S=1 

ENDIF 

IF ( S. EQ. 1 ) GO  TO  60 

IF!  (SP. EQ.  'S  '  ) . OR.  !SP. EQ.  'D '. AND. CDL . NE.  'A ' )  ) THEN 
WRITE! 1  .  14 ) 

WR  I  TE  (  5.  14  ) 

14  FORMAT! 1  OX.  *  NORMAL  TIGHTENED 

%  'REDUCED  *') 

WR 1 TE ! 1 . 16) 

WR  I  TF  !  5.  16) 

16  FORMAT! 10X,  54!  '-')  ) 

WRITE!  1.  1 0 ) NFA. TF  A, RFA 
WRITE!  5,  1  O  )  NFA,  TFA,  RFA 

10  FORMAT!  12X,  'AC  1  =  ' ,  I  3,  1  OX  ,  '  AC  1  ~  ' ,  17,  1  OX  ,  '  AC  1-  ',  17) 

WR I TE !  1 ,  20  >  F  NR , F  TR ,  FRR 


I 
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A 


t.ri  I  I 


01  7  2 
01  73 
0  1  7  'l 
017  0 


0200 
0201 
0202 
0203 
0204 
0205 
0206 
-0207 
0208 
0209 
0210 
021  1 
0212 
0213 
0214 
0215 
0216 
0217 
0218 
0219 
0220 
0221 
0222 
0223 
0224 
0225 
0226 
0227 


WR  I  IF  i  5.  20  )  I  Nl 
FORMAT  (  1 2X  .  ' Rt 
CMP  I ! 

ir  (SR.  i  o.  - D  . , 


4U;.  C  Di¬ 


01  76 

WRI  TE  (1,  33  >18.  KS  RTS 

01  77 

WR  I  7  L  <  3,  2  5  )  F  S ,  FS,  RFS 

0170 

23 

FORMAT (12X,  'Nl  15,  10X,  Ml 

8  19, 

lOX,  Nl  ' 

19) 

0179 

WRITE  ( 1, 28) 

0  1  80 

WR  I  TE  (  5,  20) 

0181 

28 

FORMAT (  10X,  34 (  '  )  ) 

0102 

END  IF 

0183 

IFtSP.EQ.  'D'.  AND.  CDL.  NE.  'A 

'  )  THEN 

0104 

WRITEC 1, 30 ) NTN, TN,  7 

0105 

WR  I  TF  (  5,  30  )  NTN,  TN,  Z 

0136 

30 

FORMAT  (  1  OX,  'AC  2=',  13,  JOX, 

'AC  2- 

,  17,  JOX. 

AC  2 

0107 

WR  I  TE  (  1,  40  )  SNR  ,  STR,  SRR 

0180 

WRITE( 5, 40) SNR, STR, SRR 

0  l  89 

40 

FORMAT  (  1  2X ,  'RE  2=',  13,  10X. 

'RL  2 

'  .  17,  10X, 

'  RE:  2 

0190 

WR  I  TE  (  1,  45  )FS,  FS,  RFS 

0191 

WR  I  TE  (5,  45)FS,  FS,  RFS 

0192 

4  5 

FORMAT!  12X,  ' N2r-  ',  15,  lOX,  'N2  ',  19 

.  lOX,  N2: 

',  19) 

0193 

END  IF 

01  94 

WRITE ( 1 , 40) 

0195 

WR  I  TE  (  5,  48V 

0196 

48 

FORMAT( 10X,  54 (  '  -  )  ) 

0197 

IF (SP. EQ.  'S'  ) THEN 

0198 

WR I TE ( 1 , 30 ) NQ, NQ, RN 

0199 

WR  I  TE  (  5,  50  )  NQ,  NQ,  RN 

ll  0<  10X,  ' N-- 


. AND. CDL. EQ.  '  A  ' .  AND.  J.  GE.  19).  OR.  (SP.EQ.  'D  ' .  AND.  CDL. 


)  FORMAT  (  1  2X .  'N~  ' ,  16,  10X.  'N=',  110,  10X,  '  N;-  ' ,  110) 

WR I TE ( 1 ,  5 1  ) 

WR I TE  <  5,  51  ) 
l  FORMAT ( 1  OX,  54  (  ' )  ) 

END  IF 

IF(  (SP.  EQ.  'D'.  AND.  CDL.  EQ.  'fit' .  AND.  J.  GE.  19).  OR.  (SP.EQ.  'D'.  A 

*  EQ.  'fit'.  AND.  M.  LE.  16)  )  THEN 
WRITE ( 1 , 54) 

WRITE  (5,  54) 

1  FORMAT ( 10X, 'SAMPLE  SIZE  NOT  AVAILABLE  AT  ANY  INSPECTION', 

*  '  LEVEL',/,  'CORRESPONDING  SINGLE  SAMPLING  PLAN  WILL', 

*  'BE  USED: ' ) 

SP=  'S' 

GO  TO  1 
END  IF 

PRINT*, 'DO  YOU  WANT  SCHEME  EVALUATION..?' 


PRINT  *,  'IF  YES  ENTER . 1' 

PRINT  *.  'IF  NO  ENTER . 2' 

READ  *, AZ 

IF(AZ. EQ. 2)  GO  TO  60 
IF( AZ. EQ. 1 )  GO  TO  58 
NNNN=L 

IF ( SP. EQ.  'S' )  JT  = 1 
IF (SP. EQ.  'D  '  )  JT -2 

CALL  OC  <  NNNN,  FS,  RFS,  RFA,  SRR ,  FRR .  NQ,  RN,  TFA,  TN,  FTR ,  STR,  NFA,  NTN 
,  FNR ,  SNR,  Z,  JT,  AQL  ) 

RETURN 

END 
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-UIJRDU  I  INI 


.)()(<  1 
ooo;> 
0003 

GOO') 
0005 
0006 
0007 
0008 
0007 
0010 
OOl  1 
0012 
0013 
0014 
0015 
OO  1  6 
0017 
0010 
0017 
0020 
0021 
0022 
0023 
0024 
0025 
0026 
0027 
0023 
0027 
0030 
0031 
0032 
0033 
0034 
0035 
0036 
0037 
0038 
0037 
0040 
004  1 
0042 
0043 
0044 
0045 
0046 
0047 
0043 
0047 
0050 
0051 
0052 
0053 
0054 
0055 
0056 
0057 


t,  a  )  i,  -I,  t>  1C  i-  i>  1:  i, 

C  *►  THIS  SUBP.OUT INI  F  INF'"",  THf-  SAMP.  F  '  :  -  OP 

C  »  >  DOTH  THL'  DUUBL.F  AND  1  HD  SINGLE  Sal'll,  i  NG  PLANE,  ■- 1 

C»*>  AT  NDRMA1  TIGHTENED  AND  REDUCED  INU'IHIOM  c  v. 

C«t  LEVELS.  «■* 

C*#»t<Ht*#***<HHHHIt****»**»#******#«-*HH*#******t*»**#*# 


CHARACTER  I  *2,  CDL  .  SR 

INTEGER  FS,  FA,  SA,  AR ,  SR,  AC,  AT,  AE,  SC,  ST,  Z,  NN,  NFS,  NFA,  NAE,  NSA 
INTEGER  NAC, NAT,  TN,  TFS,  TFA,  TAE.  TSA,  TAR,  TAC .  TAT, RN,  RFA,  RAE 
INTEGER  RAR, RAC, RAT, NAR, FRR , RFS,  FTR , FNR , SNR,  STR,  SRR 
COMMON/SCL  /L,  SP,  S,  I,  NQ,  T,  R,  CDL  ,  FS,  D,  J,  K,  AQL,  M,  FA,  SA,  AR 

*  -  SR,  AC,  AT  ,  AE,  SC,  ST,  Z,  NTN,  NFS,  NFA,  NAE,  NSA,  NAR, 

*  NAC,  NAT,  TN,  TFS,  TFA,  TAE,  TSA,  TAR-  1  AC  .  TAT,  RN,  RFS,  RFA 
T>  ,  RAE,  RAR,  RAC,  RAT,  FRR,  FTR,  FNR,  SNR.  STR.  SRR 

I E  < SF . EG.  'S'  ) THEN 

II-  (CDL.  EQ.  'A  '  )  THEN 


II  (CDL.  EQ.  'A 
NO  2 
L  ND I  F 

ir-  (CDL  .  EQ.  'B 
NQ-  3 
E  ND  I  E 

If-  (CDL.  EQ.  'C 
NQ--  5 
END  If 

If  (CDL. EQ.  'D 

NQ=G 

END  IF 

IF ( CDL. EQ.  'E 
NQ=  1  3 
END  IF 

IF (CDL. EQ.  'F 

NQ~20 

END  I F 

IF ( CDL. EQ. 'G 

NQ-32 

END  IF 

IF (CDL. EQ. 'H' 
NQ--  50 
END  IF 

IF ( CDL. EQ. 'J- 

NQ=80 

END  IF 

IF ( CDL. EQ.  'K  ' 
NQ= 1 25 
END  IF 

IF (CDL. EQ.  ’L  ' 

NQ-200 

END  IF 

IF (CDL. EQ.  'M  ' 
NQ=31 5 
END  IF 

I  F ( CDL . EQ.  'N ' 

NQ=500 

END  IF 


)  THEN 


THEN 


THEN 


)  THEN 


)  THEN 


)  THEN 


THEN 


)  THEN 


)  THEN 


)  T!  IFN 


)  THEN 


)  THEN 


oo:»;j 
0059 
0060 
an/,  t 

oo,'.; 
o<>63 
006-0 
006  5 
0066 
0067 
0068 
0069 
0070 
0071 
0072 
0073 
007-0 
007  5 
0076 
0077 
0078 
0079 
0080 
0081 
0082 
0083 
0084 
0085 
0086 
0087 
0088 
0089 
0090 
0091 
0092 
0093 
0094 
0095 
0096 
0097 
0098 
0099 
0100 
0101 
0102 
0103 
0104 
0105 
0106 
0107 
0108 
0109 
0110 
01 1  1 
0112 
0113 
0114 


ii  (cm. .  ro.  p  '  niu-N 
NG;  800 
END  IF 

If  (GDI  .  1:0.  O')  THEN 

NO--  1250 
END  IF 

IF (CDL.  EG.  'Ft  '  )  THEN 

NQ-2000 

END  I F 

IF (CDL.EQ.  'S'  ) THEN 
NO-3  1 50 
END  1 F 

IF (CDL.EQ.  'A ' . OR. CDL. EQ.  '  B  ' . OR. CDL. EQ.  'C'lTHEN 

RN-2 

END1F 

IF"  (CDL.  EQ.  'D  '  )  TFtEN 
RN~  3 
I  NDIF 

IF (CDL. EQ.  'E  '  )  I  HEN 

RN=5 

ENDIF 

IF (CDL. EQ.  'F ' ) THEN 

RN=8 

ENDIF 

IF (CDL. EQ.  'G  '  ) THEN 

RN=  1 3 

ENDIF 

IF  <  CDL.  EQ.  'H  '  ntlEN 

RN=20 

ENDIF 

IF  (CDL.  EQ.  '  J'  )  THEN 

RN=32 

ENDIF 

IF(CDL. EQ.  'K  '  ) THEN 

RN=50 

ENDIF 

IF ( CDL. EQ.  'L  '  )THEN 

RN=80 

ENDIF 

IF ( CDL. EQ.  'M '  ) THEN 

RN=125 

ENDIF 

IF(CDL. EQ.  'N  '  ) THEN 

RN=200 

ENDIF 

IF (CDL. EQ. 'P'lTHEN 

RN=31 5 

ENDIF 

IF (CDL. EQ. 'Q ' ) THEN 

RN=500 

ENDIF 

IF (CDL. EQ. 'R ' ) THEN 

RN=800 

ENDIF 

ENDIF 

IF (SP. EQ.  D'. AND. CDL. EQ. 'A' ) THEN 
FS=0 
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Ul  1.1 
0116 
0117 

o  i  i  n 

0117 
0120 
0121 
0122 
0123 
0 1  24 
0125 
0126 
0127 
0128 
0129 
0130 
0131 
0132 
0133 
0134 
0135 
0136 
0137 
0138 
0139 
0140 
0141 
0142 
0143 
0144 
0145 
0146 
0147 
0148 
0149 
0150 
0151 
0152 
0153 
0154 
0155 
0156 
0157 
0158 
0159 
0160 
0161 
0162 
0163 
0164 
0165 
0166 
0167 
0168 
0169 
0170 
0171 


NL.  A 


!  Ill  .N 


L  ND 1 1 

lMSP-EQ.  D  .  AND.  CDL  , 

If  (COL.  FQ.  '0  ITHFN 
PS  2 
ENDIf 

IP  (  CDL  .  EQ.  C  '  )  THEN 
FC  3 
END  IP 

IF (CDL. EQ.  'D  '  ) THEN 
F3—5 
END  I F 

IF ( CDL. EQ. 'E ' ) THEN 
FS=8 
END  IF 

I F ( CDL. EQ.  'F  '  I THEN 
F  S-:  1  3 
END  IF 

IF (CDL. EQ.  'G  '  ITHFN 
FS=20 
END  IF 

IF (CDL. EQ.  'H  'ITHFN 
F5~32 
END  IF 

IF (CDL. EQ.  'J'  I  THEN 
FS-  50 
END  IF 

I F ( CDL. EQ.  'K  '  ITHPN 
FS=80 
ENDIF 

IF ( CDL. EQ.  'L  'ITHFN 
FS=1 25 
ENDIF 

IF  <  CDL. EQ.  'M ' 1THEN 
FS=200 
ENDIF 

IF<CDL. EQ.  'N'lTHEN 
FS=31 5 
ENDIF 

IF (CDL. EQ.  'P ' 1  THEN 
FS=500 
ENDIF 

IF (CDL. EQ.  'Q ' 1  THEN 
FS=800 
ENDIF 

IF ( CDL. EQ.  'R ' >  THEN 
FS=1250 
ENDIF 

IF (CDL. EQ. 'S' 1THEN 
FS=2000 
ENDIF 

ENDIF 

IF ( (SP. EQ.  'D  '  I. AND.  (CDL. EQ.  'B  ' 

RFS=0 

ENDIF 

IF( SP. EQ.  'D '.  AND.  CDL. NE.  ' B ' . AND. CDL . NE.  'C  '  I  THEN 
IF (CDL. EQ. 'D' ITHFN 
RFS=2 


OR. CDL. EQ.  'C  '  1 ITHEN 


► 

5  -*• 

• 

oi  /;> 

LNim 

*  . 

t  r '  } 

If  (CDI  FQ.  T'HHf-N 

. 

0  1  7  4 

RF8  3 

1 

0  1  7  0 

ENDIF- 

im 

0176 

IF (CDL. EQ.  F  t  THEN 

. 

Ol  77 

Rh'3=5 

,v. 

0178 

END  IF 

0179 

I F ( CDL . EQ.  'G  '  ) THEN 

0180 

RFS~8 

0181 

ENDIF 

- 

0182 

I F  <  CDL.  EQ.  'H  '  )TlfEN 

J  ' 

i 

0183 

RFS= 1 3 

0184 

ENDIF 

* 

0185 

IF (CDL. EQ.  ' J  '  ) THEN 

*» 

0186 

RFS=20 

0187 

ENDIF 

0188 

I F  (  CDL.  FQ.  'K  '  )THFIN 

0189 

RFS-32 

0190 

ENDIF 

Za 

0191 

IF (CDL. EQ.  'L  '  ) THEM 

0192 

RFS=50 

0193 

ENDIF 

0194 

IF  (CDL.  EQ.  '  M  '  )THF:N 

0195 

RFS=80 

0196 

ENDIF 

' 

0197 

IF (CDL. EQ.  'N'  ) THEN 

0198 

RFS= 1 25 

m 

0199 

ENDIF 

0200 

IF ( CDL. EQ.  'P  '  ) THEN 

* 

0201 

RFS=200 

0202 

ENDIF 

0203 

IF ( CDL. EQ.  'Q  ' ) THEN 

0204 

RFS=31 5 

0205 

ENDIF 

0206 

IF ( CDL. EQ.  'R  '  ) THEN 

3 

0207 

RFS=500 

■ 

0208 

ENDIF 

- 

0209 

ENDIF 

0210 

RETURN 

021  1 

END 

a 

— 

■ 

• 

;v 

• 
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OOOl 

0002 

0003 

0004 
00 00 
0000 
0007 
0008 
0009 
OOIO 
OOl  1 

0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 

0021 

0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 


C  ***««***«♦**»**«#»*#**#»#*♦***«»»***«**##*»**#♦#«« 
C**  THIS  Si  in  ROUTINE  LOCATES  THE  ^ELl  IN  T1H  ** 
C**  f  ABLE S  PRESENTED  IN  MIL-STD-105D  CONTAIN!  ** 
C**  ACCEPTANCE  AND  REJECTION  NUMBERS  FOR  ANY  *# 
C«*  PATICULAR  COMBINATION  OF  AQL  AND  CODE  LETTER.  ** 
C*************-****************  **************** ****** 
C 


CHARACTER  1*2,  CDL,  SP 

INTEGER  FS,  FA,  SA,  AR,  SR,  AC,  AT,  AE,  SC,  ST,  l 

COMMON/SCL/L,  SH.  S,  I,  NQ,  T.  R,  CDL,  FS,  D,  J,  K.  AQL,  M,  FA.  SA,  AR 
*  SR,  AC,  AT,  AE,  SC,  ST,  Z 
IF  ( CDL  .  EQ.  'AM  THEN 
K=  1 

ELSF  IFTCDL.EQ.  'BMTMEN 
H=-2 

ELSE  IF (CDL. EQ.  C  M THEN 
K~  3 

ELSE  IFTCDL. EQ. 'D M THEN 
K-~4 

ELSE  IFICDL.  EQ.  'E  MTHEN 
K-5 

ELSE  IF(CDL.EG.  'FMTHEN 
K'_6 

ELSE  IF1CDL.EQ.  'GMTHEN 
K-  7 

ELSE  IF<CDL.EQ.  'HMTHEN 
K=8 

ELSE  IF<CDL.  EQ.  'JMTHEN 

ELSE  IFtCDL.  EQ.  'K  MTHEN 
K=10 

ELSE  IF( CDL.  EQ.  'L  MTHEN 
K=1 1 

ELSE  IF(CDL.EQ.  'MM THEN 
K—  1 2 

ELSE  IF1CDL.EQ.  'N MTHEN 
K~  1 3 

ELSE  IFICDL.EQ.  'P  MTHEN 
K=14 

ELSE  IF ( CDL.  EQ.  'Q  MTHEN 
K=  1  5 

ELSE  IF<CDL.  EQ.  'R  MTHEN 
K=16 

ELSE  IF ( CDL.  EQ.  'S  MTHEN 

K=17 

END  IF 

IF ( AQL. EQ. O. 010 ) THEN 

J-l 

ELSE  IF ( AQL. EQ. 0. 015) THEN 
J=2 

ELSE  IF ( AQL. EQ. . 025) THEN 

J-3 

ELSF  IF <AQL. EQ. 0. 040)THEN 
J=4 

ELSE  IF ( AQL. EQ. 0. 065) THEN 
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0008 
0009 
0060 
006  t 
0062 
0063 
0064 
0060 
0066 
0067 
0068 
0069 
0070 
0071 
0072 
0073 
0074 
0070 
0076 
0077 
0078 
0079 
0080 
0081 
0082 
0083 
0084 
0080 
0086 
0087 
0088 
0089 
0090 
0091 
0092 
0093 
0094 
0095 
0096 


1.1  SI  it  «  Aiii  .  EU.  0.  10;  Till  N 
J  6 

n  SF  H  'AGI  FG.  n  15)  THEN 
J— 7 

ELSE  ir  (AQL.  Ed.  0.  25)  THEN 
J-  8 

ELSE  1 F <  AQL. EG. 0. 4 ) THEN 
J=9 

ELSE  IF ( AQL. EQ. 0. 65) THEN 
J^IO 

ELSE:  IF<  AQL.  EQ.  1  )THEN 
Jr-  1  1 

ELSE:  IF(  AQL.  EQ.  1.  5)  THEN 
J-  12 

ELSE  IF  ( AQL.  EQ.  2.  5  )  THFN 
J  1  3 

ELSE  IF  (AQL.  EQ.  4)  THEN 
J  1  4 

El  SE  IF  (AQL.  EQ.  6.  5)  THEN 
1  5 

ELSE  1F< AQL. EQ. 10) THEN 
J- 1  6 

ELSE  IF(AQL.EQ. 15) THEN 
J-  1  7 

ELSE  IF  ( AQL.  EQ.  20)  THEN 
J=  1  3 

ELSE  IF(  AQL.  EQ.  40)  THEN 
J=19 

ELSE  IF ( AQL. EQ. 65) THEN 
J-20 

ELSE  IFLAQL. EQ. 100JTHEN 
J=21 
ELSE 
0=22 

END  IF 
M=J+K 
RETURN 
END 
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0001 
0002 
0003 
0004 
0005 
0006 
0007 
0008 
0009 
0010 
OOl  1 
0012 
0013 
0014 
0015 
0016 
0017 
0018 
0019 
0020 
0021 
0022 
0023 
0024 
0025 
0026 
0027 
0028 
0029 
0030 
0031 
0032 
0033 
0034 
0035 
0036 
0037 
0038 
0039 
0040 
0041 
0042 
0043 
0044 
0045 
0046 
0047 
0048 
0049 
0050 
0051 
0052 
0053 
0054 
0055 
0056 
0057 


SUBROUTINE  V/M  utlfj 
C 

C**  »*«»l«***t**»*««««*******#«**«*»*«#*»#**ft*«»««*«*IH>  to  (Mm** 

C**  THIS  SUBROUTINE  F INDS  ACCFPTANCF  AND  REJECTION  »# 

C**  NUMBERS  FOR  THE  CELL  IDENTIFIED  BY  SUBROUTINE  INDEX.  ** 

C********  **#•**#■»•*  **-**#**#******#***»********-tt********#*-l<-*-M-*« 

c 

CHARACTER  1 *2. CDL. SP 

INTEGER  FS, FA,  SA, AR ,  SR, AC, AT, AE, SC, ST,  Z.  NN,  NFS.  NFA.  NAE,  NSA 
INTEGER  NAR, NAC. NAT,  TN. TFS, TF A,  TAE. TAR, TAC,  TAT.  RN,  RFA,  RAE 
INTEGER  FRR,  FTR , FNR , SNR, STR,  SRR,  RAR,  RAC,  RAT,  RFS,  TSA 
CQMMON/SCL/L.  SP,  S,  I,  NQ,  T,  R,  CDL,  FS,  D,  J,  K,  AQL,  M.  FA,  SA. 

♦  AR,  SR,  AC,  AT,  AE,  SC,  ST,  Z,  NTN,  NFS,  NFA,  NAE,  NSA,  NAR.  NAC,  NAT,  TN 
S  .  TFS,  TFA.  TAE,  TSA-  TAR,  TAC.  TAT,  RN,  RFS,  RFA-  RAE,  RAR,  RAC.  RAT 

*  ,  FRR, FTR,  FNK, SNR,  STR,  SRR 
Nl  A—  1  00 

ti-'a=  i  oo 

RFA= 1 00 
NAT  = 1 00 
TAT=100 
RAT=100 
TAC  =3  00 
RAC= 1 OO 
T  AR- 1 00 
RAR= 1 00 
T  AE- 1 00 
RAE- 1 00 
TSA=100 
NAC=  1 00 
NAR-'  1  OO 
NAE= 1 00 
NSA— lOO 
STrrlOO 

sc=ioo 

SR- 100 

IFISP.EQ.  ‘S ‘ )THEN 

IF  <  M. GE. 2. AND. M. LE.  16>NFA=0 
IF(M. GE. 26)NFA=2l 
IF  <  M. EQ. 25>NFA=14 
IF ( M. EG. 24 ) NFA= 10 
IF ( M. EQ. 23 ) NFA=7 
IF(M. EQ. 22)NFA=5 
IF(M. EQ. 21 ) NFA=3 
IF ( M. EQ. 20 ) NFA— 2 
IF ( M. EQ.  19 ) NFA= 1 
IF  (  M.  EQ.  17.  AND.  J.NE.  16>NFA=0 
IF  <  M. EQ.  18. AND. J. NE.  18)NFA=1 
IF ( K. EQ.  16. AND. J. LE.  2)NFA=0 
IF( J. EQ. 16. AND. K. LE. 2)NFA=1 
I F  <  M. GE. 26 ) TFA= 1 8 
I F ( M. EQ. 25 ) TFA=12 
IF ( M. EQ. 24 ) TFA=8 
I F ( M. EQ. 23 ) TFA=5 
IF  ( M.  F.Q.  22  )  TFA=3 
IF (M. EQ. 21 )TFA-2 

IF  ( M.  GE.  18.  AND.M.LE.  20.  AND.  J.NE.  16>TFA=1 
IF (K. EQ. 16. AND. J. EQ. 2 ) TFA=2 


0050 
oo  :>v 
0060 
006  J 
0062 
0063 
0064 
0065 
0066 
0067 
0068 
0069 
00/0 
0071 
0072 
0073 
0074 
0075 
0076 
0077 
0070 
0079 
0080 
0081 
0082 
0083 
0084 
0085 
0086 
0087 
0088 
0089 
0090 
0091 
0092 
0093 
0094 
0095 
0096 
0097 
0098 
0099 
OlOO 
OlOl 
0102 
0103 
0104 
0105 
0106 
0107 
0108 
0109 
Ol  lO 
Olll 
0112 
0113 
0114 


Jt  (I'l.  I  E.  1  7.  AND.  J.  LI.  .15)  1  f  A  O 
IF  ( J.  HO.  16.  ANL).  K.  IJ.  2  )  IT  A  I 
IF  <1*1. GE. 26) RAC- 10 
IF  <1*1. EG. 25 ) RAC- 7 
IF (M. HQ. 24 ) RAC— 5 
IF (M. HQ. 23  >  RAC-3 
IF (M. HQ. 22. AND. J. LE. 191RAC -2 
IF  <  1*1. HQ. 21 ) RAC  =  1 
IF (M. EQ. 20. AND. J. LE. 18)RAT=1 
IF  <1*1. EQ. 19. AND. J. LE. 17>RAT^0 
I F  ( 1*1.  GE.  2.  AND.  1*1.  LE.  I6)RFA=0 
IF  ( 1*1 .  EQ.  17.  AND.  J.  LE.  I5)RFA=0 
IF  < J. LE. 2)RFA=0 
IF ( J. EQ. 16. AND. K. LE. 3>RAT=0 
IF  ( J.  EQ.  17.  AND.  K.  EQ.  1  >RAT=0 
IF  < J. EQ. 21. AND. K. LE. 2)RFA=5 
IF ( J. EQ. 20. AND. K. EQ. 1 >RFA=3 
IF ( J.  EG.  20. AND. K. EQ. 2)RAT=3 
IF  <  J.  EQ.  19.  AND.  K.  EQ.  1  )RFA=2 
IF  <  J. EQ. 19. AND. K. EQ. 2)RAT=2 
IF ( J. EQ. 18. AND. K. EQ. 1 >RFA=1 
IF ( J. EQ. 18. AND. K. EQ. 2)RAT=1 

ENDIF 

IF(SP. EQ.  'D' ) THEN 

I F ( M. LE. 16. AND. J. LE. 15>THEN 

NAT-0 

NSA- 1 

TAT=0 

TSA- 1 

RAT=0 

ST=0 

END  IF 

IF  <K. EQ. 16. AND. J. LE.2)THEN 
NAT=0 
NSA=1 
TAT=0 
TSA=1 
RAT =0 
ST =0 
ENDIF 

IF  <  K.  EQ.  1.  AND. J. EQ.  18) THEN 
NSA=3 
NAC=0 
TAT=0 
TSA=  1 
ENDIF 

IF  <M. GE.  26) THEN 

NAE=1 1 

NSA=26 

ENDIF 

IF  < M. EQ.  25 ) THEN 
NAR~7 
NSA=18 
ENDIF 

IF  <  M. EQ. 24) THEN 
NAR-5 
NSA= 1 2 
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VAI  l  f|  . 

Oil  5 

iNim- 

0  1  1  6 

I F ( M. EQ. 23)  f  HEN 

Oil/ 

NAR-3 

O'  1!' 

NSA—8 

0  11  V 

END1F 

0120 

I  F  ( M. EG. 22) THEN 

0121 

NAC-2 

0122 

NSA-6 

0123 

END  IF 

0124 

IF (M. EG. 21 ) THEN 

0125 

NAC-l 

0126 

N8A=4 

0127 

TNDIF 

0128 

I F  <  M. EQ. 20) THEN 

0127 

NAC=0 

0130 

NSA=3 

0131 

ENDIF 

0132 

If  ( (M. EQ.  17.  AND. J. GT. 2) . OR.  (M. EQ.  18. AND. J. GT. 2) ) THEN 

0133 

NAT-0 

0134 

NSA- 1 

0135 

ENDIF 

0136 

I F ( J. EQ.  16.  AND.  K.  LF. 3 ) T  HEN 

0137 

NAT=0 

0138 

NSA^l 

0 1  37 

ENDIF 

0140 

IF  ( M. GE. 26) THEN 

0141 

TAE=7 

0142 

TSA=23 

0143 

ENDIF 

0144 

IF(M. EQ. 25) THEN 

0145 

TAR  =6 

0146 

TSA=1 5 

0147 

ENDIF 

0148 

IF(M. EQ.  24) THEN 

0147 

TAR=3 

0150 

TSA~ 1 1 

0151 

ENDIF 

0152 

IF(M. EQ.  23) THEN 

0153 

TAC=2 

0154 

TSA=6 

0155 

ENDIF 

0156 

IF (M. EQ.  22) THEN 

0157 

TAC=1 

0158 

TSA=4 

0157 

ENDIF 

0160 

IF ( M. EQ. 21 ) THEN 

0161 

TAC=0 

0162 

TSA=3 

0163 

ENDIF 

0164 

IT  <M.  GE.  18.  AND.  M.  LE.  20.  AND.  J.  GE.  3)  THEN 

0165 

TAT=0 

0166 

TSA=  1 

0167 

ENDIF 

0168 

IF(M. EQ. 26 ) THEN 

0167 

RAE  =  5 

0170 

SR=  1  2 

0171 

END  If 

’  Ai  >  •;  ' 

0 ;  72 

II  <H.  F Q .  25)  1  HI  N 

01  73 

waf  :t 

0174 

Sl<~0 

0175 

r-Noif 

0!  76 

I F  ( M .  EQ.  24  )  THEN 

0177 

RAK=2 

0170 

SC =6 

0177 

END  11 

01  BO 

IF  ( M. EQ. 23) THEN 

0181 

R  AR~- 1 

0182 

SC=4 

0183 

END  IF 

0184 

I F ( M. EQ. 22  >  THEN 

0185 

RAR-0 

0186 

SC-3 

0187 

END  IF 

0188 

I F  <  M. EQ. 21 ) THEN 

0189 

RAR=0 

0190 

SR  =  1 

0191 

ENDIF 

0192 

IF  <  M. EQ. 20  >  THEN 

0193 

RAC=0 

0194 

SR~0 

0195 

ENDIF 

0196 

IF  (M.  EQ.  1  9  )  THEN 

0197 

RAT=0 

0198 

ST -0 

0199 

ENDIF 

0200 

IF  (  M.  EQ.  18.  AND.  J.  GE 

0201 

RAT=0 

0202 

ST=0 

0203 

ENDIF 

0204 

ENDIF 

0205 

RETURN 

0206 

END 

15) THEN 
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hubhoutini  cum 


OOO  i 
0002 
0003 
0004 
0005 
0006 
0007 
0008 
0009 
0010 
OOll 
0012 
0013 
0014 
0015 
0016 
0017 
0018 
0019 
0020 
0021 
0022 
0023 
0024 
0025 
0026 
0027 
0028 
0029 
0030 
0031 
0032 
0033 
0034 
0035 
0036 
0037 
0038 
0039 
0040 
0041 
0042 
0043 
D044 
0045 
D046 
0047 
0048 
0049 
0050 
0051 
0052 
0053 
0054 
0055 
0056 
0057 


C 

C »***##»«*#*##«  «***«*t»***»«'«#**«'««»*««*t*tt#ft«#*tt**ii*«« •***•**#* 

C**  TH I P  SUBROUTINE  IDENTIFIES  THE  CODE  LETTER  ASSOCIATED  ** 
C**  WITH  EITHER  THE  SPECIAL  OR  GENERAL  INSPECI10N  LEvLlS  **• 
C*«  FOR  ANY  GIVEN  LOT  OR  BATCH  SIZE.  ** 

C  ******************************************************  *■**•***•* 

c 


CHARACTER  1*2. CDL. SP 
INTEGER  S1.S2,  S3,  S4 
COMMON/SCL /L.  SP.  S,  I ,  NO,  T.  R,  CDL 
IF ( I. EQ.  'SI  '  ) THEN 

IF (L. GE. 2. AND. L. LE. 50) CDL ^  'A ' 
IF(L. GE.  51. AND. L. LE. 500)CDL='B ' 

IF ( L.  GE.  501. AND. L. LE. 35000 ) CDL  =  'C 
IF(L. GE. 35001 ) CDL= 'D ' 

END  IF 


IF ( I. EQ.  ' S2 ' ) 7  HEN 


IF(L.  GE.  2.  AND.  L  .  LE.  25)CDL-^  'A' 

IF ( L. GE. 26. AND. L. LE. 1 50 ) CDL= ' B ' 

IF  <L. GE.  151. AND.L.LE.  1200>CDL='C' 

IF ( L. GE. 1201. AND. L. LE. 35000 ) CDL= ' D ' 

IF  <  L. GE. 35001 )CDL='E' 

END  IF 

I F ( I. tO.  'S3'  ) THEN 

IF  <  L. GE. 2. AND. L. LE.  15)CDL=  'A' 

IF (L. GE. 26. AND. L. LE. 50)CDL=  'B  ' 

IF ( L. GE. 51. AND. L. LE. 150)CDL= 'C ' 

IF ( L. GE. 151. AND. L. LE. 500,CDL='D' 

IF vL.  GE.  501 . AND . L . LE. 3200 ) CDL= ' E ' 

IF (L. GE.  3201 . AND. L. LE. 35000 )CDL= 'F ' 
IF(L. GE. 35001 . AND. L. LE. 500000) CDL = 'G ' 
IF(L. GE.  500001 >CDL='H' 

END  IF 

IF(I.EQ.  'S4  '  )  THEN 

IF <L.  GE.  2. AND. L. LE. 15)CDL=/A' 

IF (L.  GE.  16.  AND.  L. LE. 25>CDL='B ' 

IF <L.  GE. 26.  AND. L. LE. 90)CDL=  'C  ' 
IF1L.GE.  91.  AND.L.LE.  150>CDL='D' 

IF (L.  GE.  151.  AND. L. LE. 500)CDL='E' 

IF(L. GE. 501. AND. L. LE. 1200)CDL='F' 

IF  < L.  GE.  1201. AND. L.LE.  10000)CDL=  'G ' 

IF ( L.  GE.  10001.  AND.L.LE. 35000 ) CDL= 'H ' 
IF ( L.  GE.  35001.  AND. L. LE. 500000 ) CDL= 'J ' 
IF ( L.  GT.  500000 ) CDL—  'K  ' 

END  IF 

IF ( I. EG. '1 ' ) THEN 

IF(L. GE. 2. AND. L.  LE.  15)CDL=,A' 

IF ( L. GE. 16. AND. L. LE. 25)CDL= 'B ' 

1 F ( L . GE. 26. AND. L. LE. 90)CDL='C ' 

IF(L. GE. 91. AND.  L.  LE.  150)CDL= 'D' 

IF (L. GE. 1 51. AND. L. LE. 280 )CDL= 'E ' 

IF (L. GE. 281. AND. L. LE. 500)CDL= 'F ' 
IF1L.GE. 501. AND. L. LE. 1200)CDL=/G' 
IF(L. GE. 1201. AND. L. LE. 3200)CDL='H' 

IF ( L. GE. 3201. AND. L. LE. 1 0000 ) CDL= ' J ' 

IF ( L. GE. 10001. AND. L. LE. 35000 )CDL= 'K ' 


oodh 
oo  59 
0060 
0061 
006? 
0063 
0064 
0065 
0066 
0067 
0068 
0069 
0070 
0071 
0072 
0073 
0074 
0075 
0076 
0077 
0078 
0079 
0080 
0081 
0082 
0083 
0084 
0085 
0086 
0087 
0088 
0089 
0090 
0091 
0092 
0093 
0094 
0095 
0096 
0097 


It  (l_.  Gfi  .  35001.  AND.  L.  LE.  1  50000  )  CDL-  '  L  ' 

It  <L. GE.  150001. AND. L. LE.  500000) CDL- ' M ' 

IF (L. GE. 500001 )CDL- 'N ' 

LNP!  F 

II  (  I  .  t*G.  '2  '  )  THEN 

IF (L.  GE.  2. AND. L . LE. 3) CDL -  'A  ' 

IF  <L. GE. 9. AND. L. LE.  1 5 ) CDL-  ' D  ' 

IF(L. GE. 16. AND. L. LE. 25) CDL- 'C ' 

IF (L. GE. 26. AND. L.  LE.  50)CDL= 'D ' 

IF ( L. GE. 51. AND. L. LE. 90JCDL-  'E  ' 

IF  <  L. GE. 91. AND. L. LE.  1 50  >CDL=  'F ' 

IF ( L. GE. 151. AND. L. LE. 280 ) CDL- 'G' 
IF<L. GE. 281. AND. L . LE.  500>CDL=  'H' 
IF(L. GE. 501. AND. L. LE. 1200) CDL- ' J' 

I F ( L. GE.  1201. AND. L. LE. 3200)CDL='K  ' 
IF(L.  GE. 3201 . AND. L. LE.  10000 ) CDL=  'L 
IF(L. GE.  10001. AND. L. LE. 35000 )CDL=  '  M 
IF  (  L.  GE.  35001.  AND.  L.  LE.  1  50000 )  CDL.=  ' N 
I F ( L . GE. 150001. AND. L. LE. 500000)CDL= ' 
I F ( L. GE. 500001 )CDL='Q' 

END  I F 

IF ( I. EQ.  '3  '  ) THEN 

IF-  (L.  GE.  2.  AND.  L.  LE.  B)CDL=  'B  ' 

I F ( L . GE. 9. AND. L. LE. 15)CDL='C ' 

IF  <L.  GE.  16.  AND.L.LE.25)CDL='D' 

IF ( L. GE. 26. AND. L. LE. 50JCDL- 'E ' 

IF ( L. GE. 51 . AND. L. LE. 90)CDL= 'F ' 

IF (L. GE. 91. AND. L. LE. 150)CDL='G' 

IF ( L. GE. 151. AND. L. LE. 280)CDL='H' 

IF ( L. GE. 281 . AND. L. LE.  500 ) CDL= '  J ' 

IF<L. GE. 501. AND. L. LE. 1200)CDL='K' 

IF ( L. GE. 1201 . AND. L. LE. 3200)CDL= 'L ' 

I F ( L. GE. 320 1 . AND. L. LE. 1 0000 ) CDL= ' M ' 

IF  <  L. GE.  10001. AND. L. LE. 35000)CDL='N ' 
IF ( L. GE. 35001. AND. L. LE.  1 50000 )CDL=  'P  ' 
IF< L. GE.  150001. AND. L.  LE.  500000 ) CDL= ' Q ' 
IF  <  L.  GE. 500001 >CDL='R' 

END  IF 

RETURN 

END 
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OOO  l 

0002 

ooo:j 

0004 

oooo 
0006 
0007 
0008 
0009 
OOIO 
OOl  1 
0012 
0013 
0014 
0015 
0016 
0017 
0018 
0019 
0020 
0021 
0022 
0023 
0024 
0025 
0026 
0027 
0028 
0029 
0030 
0031 
0032 
0033 
0034 
0035 
0036 
0037 
0038 
0039 
0040 
0041 
0042 
0043 
0044 
0045 
0046 
0047 
0048 
0049 
0050 
0051 
0052 
0053 
0054 
0055 
0056 
0057 


SUliKuu  I  1 NE  OCCNNNN.  ITS..  IRFS-  IRPA.  I3RR,  IFRR.NO. 
1  IN,  If  1  R.  I STR ,  NFA,  NTN,  IFNR,  ISNR,  1  Z,  J,  AQl.  ) 


1RN,  1 TF  A 


f;*#«»*»##*****#*###*##jHHt»***#t**»***#***#***HHI***»t#»#**»**** 

C*»  IHlb  SUBROUTINE  EVALUATES  THE  SCHEME  OPERATING 

C*«  CHARACTER  1ST  ICS, FOR  THE  REQUIRED  SAMPLING  PLANS.  #* 

c 

REAL*8  AQ(20,  20),  BQ(20,  20) 

DIMENSION  N ( 30 ) # PR ( 30 ) , V< 100,  50), P (30) ,  TPA<30,  4).  ASN(30,  4) 
DIMENSION  AOQ ( 30,  4), ATI (30, 4), NN<4>, TM<20,  20),  ZP(20, 20) 
DIMENSION  ASNC (30) , AOQC (30), AT I C ( 30 ) ,  AFI ( 30 ) ,  PA ( 30 ) 

INTEGER  A(50).  R<50),  HI,  A  A  (2,  4),  RR(2,  4) 

IF ( J. EQ. 2 )  GO  TO  197 
IF  (  U.  EQ.  1  )  GO  TO  172 
1/2  NN ( 1 ) =NQ 
NN ( 2 ) =NG 
NN ( 3 ) — IRN 
NN ( 4 ) — I RN 
AA ( 1, 4 ) - I FRR- 1 
RR ( 1 , 3) = IRFA+1 
GO  TO  272 
197  NN ( 1 )=IFS 

NN ( 2 )  =  1 FS 
NN ( 3 ) =  IRFS 
NN ( 4 ) - I RFS 
AA ( 1, 4 ) = IRFA 
A A ( 2, 4 ) - ISRR- 1 
RR( 1. 3 ) =  IFRR 
272  AA(1,1)=ITFA 

AA ( 2, 1 )~ITN 
RR<  1,  1 )=IFTR 
RR ( 2, 1 )=ISTR 
AA( 1, 2)=NFA 
AA(2,  2) —NTN 
RR( 1. 2)~IFNR 
RR ( 2, 2)=ISNR 
AA ( 1 , 3>=IRFA 
AA ( 2, 3>  =  IZ 
RR (2, 3)=IZ+1 
RR ( 1 » 4) =IFRR 
RR (2, 4 )  =  I SRR 

PRINT  *,  'DO  YOU  WANT  A  TABLE  OR  A  GRAPH  FORMAT  ?' 

PRINT  *.  '  FOR  GRAPH _ ENTER:  1  ' 

PRINT  *,  '  FOR  TABLE _ ENTER:  2  ' 

READ  *,  XTC 

IF ( XTC. EQ. 1 >  GO  TO  733 

PRINT  *, 'SPECIFY  THE  NUMBER  OF  FRACTION  DEFECTIVE  VALUES.. 
READ  *,  J1 

PRINT  *,  'ENTFR  THE  FRACTION  DEFECTIVE  VALUE(S), ' 

PRINT  *,  '(PUT  A  COMMA  BETWEEN  VALUES.  ) . ' 

READ  *,  (P  (  I  ) ,  1  =  1,  J1  ) 

WRITE! 1, 27) 

WR I TE ( 5.  27) 

GO  TO  447 
733  Jl=21 

DO  55  1=1, U1 


VALUES. 


I 


O'  >58 
0059 
0060 
006  t 
0060 
006 J 
0064 
0065 
0066 
0067 
0068 
0069 
0070 
0071 
0072 
0073 
0074 
007  5 
00/6 
0077 
0078 
0079 
0080 
0081 
0082 
0083 
0084 
0085 
0086 
0087 
0088 
0089 
0090 
0091 
0092 
0093 
0094 
0095 
0096 
0097 
0098 
0099 
OlOO 
OlOl 
0102 
0103 
0104 
0105 
0106 
0107 
0108 
0109 
OllO 
011  1 
0112 
0113 
0114 


1 '  v  i  '  V  i  t  >  /',!().  O 

v.  p  (  :  ;  -  r  i  ]  '  /  if>.  o 
ao  no  t  on  i  i .  ,n 

DO  to  !  M  t  A 
01 ,  17  J  1  J 
A  l  I  AA  (  1  ,  l  M  ) 

R  (  I  >-=NR  <  I  .  1.  M  > 

N (  I  ) -NN ( I  M ) 

1/  CONTINUE 
NNN=NNNN 
M=  I  NT  <  NNN*P  ( l  )  ) 

ZNK=NNN*P  <  L  ) 

ZN=M 

IFtZNK  ZN.I.E.O.  5)  GO  TO  148 
M-M+l 

148  CONTINUE 
3UN-0 
I  1  -  R  (  1  ) 

C 

C:  .  . 

C  THE  NEXT  STATEMENTS  CALCULATE  THE  PROBABILITY  THAT 
C  THERE  ARE  (II)  DEFECTIVES  IN  THE  FIRST  SAMPLE. 

C  :  ■  :  :  : 

C 

DO  20  J -  J  »  11 

V(  I,  1) - PP ( N  < 1 ) . P  <  L ) »  I  D 
20  8UM=SUM+V< I, 1) 

V ( R ( 1 ) + 1 . 1  >  =  1. -SUM 
IF ( A < 1 ). LE. 0)  GO  TO  21 
I1=A( 1 ) 

C 

C :  :  :  ::::::::  : 

C  THE  FOLLOWING  STATMENTS  CALCULATE  THE  PROB.  THAT  THERE  : 

C  ARE  LESS  THAN  OR  EQUAL  TO  AID  DEFECTIVE  ITEMS. 

C  I:::::::;::::::::::;:::::::::;::::::::::::::::::::  t  :::::  : 

C 

DO  30  1  =  1,  II 
C 

C :::::::  : 

C  THE  NEXT  SECTION  CALCULATES  THE  INTERMEDIATE  PROBABILITIES 
C  OF  CONTINUED  SAMPLING  FOR  THE  DOUBLE  SAMPLING  PLANS. 

C  I:::::::::::::::::::::;:::::::-:::::::::::::::::::::::::::  : 

C 


30  V  <  A  <  D  + 1 ,  1) =V ( A ( 1)  + 1 ,  D  +V ( I ,  D 
21  IF ( J.  EQ.  1)  GO  TO  41 
LOW  =1 
1  =  2 

NNN=NNN-N  < 1-1) 

NR-R ( J  ) 

LOW=MAXO  <  LOW, A< 1-1) +2) 

HI=R  <  I-D 
DO  60  1 1  =  1 .  NR 
VIII, I >=0 

IF( 11. LT. LOW)  GO  TO  60 
IF(A( I ). EQ. (-1) )  GO  TO  61 
IF<A<  I  )-M.  NE.  1 1)  GO  TO  61 
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r— - -1 - - 

m  ■  «  -  4  ■  It  ■  1  ■  «  ■  JW  1  ■  J  rn  ’W  m  W  w  V  r.  V  "  ■  ■  ■  *  - - — - -  yw  . - ~ - — 

■ 

■  • 

O  1  1  ', 

MMM-M1NCH  I  1  ,  HI  :■ 

0  1  1 6 

DO  70  12-  LOW  .  MMf 

fi 

011/ 

14-11  124 1 

' 

F* 

01  IB 

DO  50  11-1, NR 

01  1  7 

PR <  I  I  ) ~PP (N< I >  r  l  ...11  ;  ■  NNN,  M  1P+1) 

0  1  20 

50 

CONTINUE 

0121 

DO  80  13=1, 14 

V- 

0  1  2P 

80 

V<  I  1,  I  )  -V  (11.  I  )  -*  V  1  T  2 ,  I  1  )  *PR  (  I3> 

jfl-  ,v 

0123 

70 

CONTINUE 

0  1  24 

GO  TO  60 

0120 

6) 

MMM=MIN0(HI,  1 1 ) 

0126 

DO  90  I2=L0W  ,  MMM 

0127 

13=11-12+1 

0120 

DO  743  IJ=1, NR 

0129 

743 

PR ( I J ) —PP ( N  < I ) ,  P ( L ) ,  IJ-1  ) 

0130 

90 

V<  I  1,  I  )=V< I 1,  I >  +  V< 12,  1-1 >*PR (13) 

0131 

60 

CONTINUE 

0132 

V(R(  I  )  +  l.  I  )  =0 

z* 

0133 

DO  110  1 1 —LOW  ,  HI 

M 

0134 

SUM=0 

0135 

I3=R (  I  )- 1 1  +  1 

0136 

DO  120  12=1, 13 

0137 

120 

SUM=SUM+PR< 12) 

0138 

1  10 

V ( R (  I  )  +  1 ,  I ) =9 ( R ( I ) + 1 ,  I )  +  ( 1. -SUM ) *V ( 1 1 ,  1-1 ) 

0139 

4  1 

CONTINUE 

0140 

C 

0141 

C 

1 

• 

0142 

C  THE 

NEXT  SECTION  EVALUATES  THE  SCHEME  FOR  THE  SINGLE 

1 

0143 

C  SAMPLING  PLAN  (DR  THE  FIRST  STAGE  IN  DOUBLE  SAMPLING). 

0144 

C 

0145 

c 

* 

0146 

121 

ASN ( L, LM ) =0 

0147 

TPA ( L, LM ) =0 

0143 

SS=0 

3 

0149 

DO  130  1  =  1,  J 

* 

0150 

SS=SS+N (  I  ) 

-  - 

0151 

IF ( A ( I ) . EG. ( — 1 ) )  GO  TO  131 

' 

0152 

ASN (L, LM ) =ASN (L,  LM )  +  (  V  ( A  (  I ) +1 ,  I )+V(R( I >  +  l,  I ) >*SS 

0153 

TPA ( L,  LM)=TPA(L,  LM > +V ( A ( I ) + 1 ,  I ) 

•.  * 

0154 

GO  TO  130 

„ 

0155 

131 

ASN ( L* LM ) =ASN ( L,  LM ) +V ( R ( I ) +1 ,  I ) *SS 

0156 

130 

CONTINUE 

3 

0157 

xxx=o 

- 

■ 

0158 

sss=o 

- 

0159 

vw=o 

0160 

DO  133  1  =  1,  J 

0161 

SSS=SSS+N ( I ) 

0162 

XXX=XXX+SSS*V< A( I >  +  l,  I) 

0163 

WV=VVV+  ( NNNN-SSS  )  *V  ( A  <  I  )  + 1 ,  I  ) 

0164 

133 

CONTINUE 

0165 

ATI  (L,  LM)=XXX  +  (NNNN4*-(  1. 0-TPA(L,  LM)  )  > 

. 

0166 

AOQ ( L, LM)  =  <  VW*P  <L  > ) /NNNN 

0167 

140 

CONTINUE 

0168 

10 

CONTINUE 

0169 

19 

FORMAT  (  10X,  F5.  3,  2X,  F9.  4,  F10.  2,  4X,  F6.  4,  IX,  FI 0.  2) 

« 

0170 

27 

FORMAT ( ///,  16X,  'SCHEME  OPERATING  CHARACTERISTICS  ',  /,  10X 

-  -  « 

0171 

* 

,'P  P ( A)  ASN  AOO  AFI') 

■-'■j 

• 

■ 

• 

r. 

r. 

• 
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.1 

1BSBI 

■  V  1 

01  77 
oi7:) 
0174 
01  75 
0  1  /A 
0177 
0170 
0179 
0180 
0181 
0181? 
0183 
0184 
0185 
0186 
0187 
0138 
01  89 
0190 
0191 
0192 
0193 
0194 
0195 
0196 
0197 
0198 
0199 
0200 
0201 
0202 
0203 
0204 
0205 
0206 
0207 
0208 
0209 
0210 
0211 
0212 
0213 
0214 
0215 
0216 
0217 
0218 
0219 
0220 
0221 
0222 
0223 
0224 
0225 
0226 
0227 
0228 


Dii  :  l  ,  70 

1)8  1  JJ  1,80 
TM< l  JJ )  O 


THIS  SFC1I0N  FIRS!  CALCULATES  THE  PROD.  Ql  DEING 
IN  NORMAL. TIGHTENED  OR  REDUCED  INSPECTION  LEVELS  AND 
THEN  COMPUTES  THE  PROG.  OF  PASSING  THE  i  IMIT  NUMBFR 
CRITERIA. 


T --  TPA  ( l  .  1  ) 

S-TPA(L.  2) 

RR 1 =TPA ( L.  3) 

R 1  -  TPA  <  L  ,  4 ) 

NZ  ~  1  0*-ASN  <  L ,  2) 

CALL  LNMIAOL,  NZ ,  L.NC  ) 
(F-IKNZ,  P  (L  )  ,  LNC  ) 


THIS  SECTION  READS  IN  THE  PROBABILITY  TRANSMION  MATRIX. 


TM  <  1 ,  1  )  =•  1 .  — T 
TM < 1 ,  2 ) -  7 
TM  (  2i  1  )■  1.  -T 
TM  (  2,  3  >  rT 
TM  <  3.  1  >  -  1 .  -T 
TM  (  3,  4  )  :-T 
TM  {  4i  1  )  -  1 .  — T 
TM  (  4,  5  )  =T 
TM ( 5.  1 )^1. -T 
TM<  5.  6  )  ~T 
TM < 6,  7)  =  1. -S 
TM ( 6, 1 1 >=S 
TM (  7.  1  >=1.  ~S 
TM  (7,8)  =S 
TM ( 8,  1  )'l.  -S 
TM  (8,9)  ~S 
TM  (  9,  1 )~1. -S 
TM ( 9,  10)=S 
TM  (  10,  1  >  =  1.  -S 
TM  < 10,  1 4 ) =S 
TM  < 1 1 , 7 )  =  1 .  -S 
TM  <11,  1  2 )  =S 
TM ( 12,  7 )  =  1 .  -S 
TM(  12,  1 3  )  =S 
TM  (  13,  7>  =  1.  -S 
TM( 13,  1 4 ) ~S 
TM ( 14,  7 ) -  1 .  S 
TM(  14,  15)  --S 
TM(  15,  7)  =  1.  -S 
TM(  15-  1 6  )  ~S 
TM  (  16,  7  )  =  1 .  -S 
TM(  16,  1  7  )  --S 
TM<  17,  7  )  —  1 .  -S 
TM  < 17,  18) =S 


0230 

r m  (  i r  ro  3 

0231 

T  M  (  1 V ,  /  )  )  . 

0232 

f  M (  1 9  19)  ( 

0333 

nr  m  <7.  30  >  o 

0234 

IN (20, 6,  1. 

0235 

TM(2C), 30)- R 

0236 

NR  -  20 

0237 

NC=20 

0238 

DO  100  1=1. 

0239 

DO  100  JJ=1 

0240 

IF ( JJ-1 ) 1 1 1 

0241 

1  5 

ZP  (  JJ,  I  >=TM 

0242 

GO  TO  100 

0243 

1  1  1 

ZP( JJ,  I ) ~  T M 

0244 

100 

CONTI NUF 

0245 

DO  115  1=2, 

0246 

DO  115  JJ-  1 

024  7 

KJ  =1  —  1 

024Q 

AQ< I,  JJ)=ZP 

0249 

1  1  5 

CONTINUE 

0250 

DO  1 8  J  J= 1 , 

0251 

AQ( 1, JU)=1. 

0252 

18 

CONTINUE 

0253 

CALL  INVERT 

02' 54 

PA ( L ) =0 

0255 

ASNC  <L )=0 

0256 

AOQC ( L ) -  0 

0257 

ATIC (L  >=  0 

0258 

DO  161  1=1,1 

0259 

I F  (  I .  GT.  5  )  Gl 

0260 

ZZ  =  T 

0261 

Z1=ASN(L,  1 ) 

0262 

Z2=AQQ ( L,  1 ) 

0263 

Z3=ATI <L> 1 ) 

0264 

GO  TO  201 

0265 

200 

IF  <  I.  GT.  19) 

0266 

ZZ=S 

0267 

Z 1 =ASN <  L,  2) 

0268 

Z2=A0Q<L,  2) 

0269 

Z3=AT I (L,  2) 

0270 

GO  TO  201 

0271 

210 

ZZ=R  1 

0272 

0273 

0274 

0275 

0276 

0277 

0278 

0279 

0280 

0281 

0282 

0283 

0284 

0285 


TO  200 


GO  TO  210 


Z1=ASN(L.  4 ) 
Z2=A0Q<L,  4) 
Z3=ATI (L,  4) 


THE  FOLLOWING  STATEMENTS  MULTIPLY  THE  PROBABILITIES 
OF  BEING  IN  THE  VARIOUS  STATES  BY  THE  CHARACTERISTICS 
OF  EACH  STATE. 


PA  (  L  )  :-PA  ( L  )  <BQ(  1,1)  *ZZ 
ASNC ( L ) =ASNC ( L  >  +BQ  < 1 ,  1 ) *Z 1 
AOQC ( L ) =AOQC ( L ) +BQ ( 1 ,  1 ) *Z2 
ATIC <L >=AT1C (L )+BQ( I.  1)*Z3 


0206 

0287 

0200 

0289 

0290 
o;?9 1 

029-7 

0292 

0294 

0295 

0296 

029/ 

0290 

0299 

0300 

0301 

0302 

0303 

0304 

0305 

0306 

0307 

0308 

0309 

0310 

0311 

0312 

0313 

0314 

0315 

0316 

0317 

0318 


A1  I  (  L.  )  A  IIC  <1  >  /NNNN 
16  1  Cl'JNT  INUt 

IF(XTC.NE.2)  GO  TO  10! 

WR  TTK  (  1  .  1  9  )  P  (  L  )  .  PA  <  l  «  ASM!.  (  L  .  .  AOUC  (  L  )  ,  A(  1(1  ) 
WRITE  (5,  19)  PH  ),  PAG  'I  i  A(J(F.  <  i  >  .  AP  I  .  I  } 

'OH  CONTINUE 

IF( XTC. EG. 2)  GO  TO  556 
WRITE! 1, 45) 

WR  I TE  (  5/  45) 

WRITE! 1, 46) 

WRITE!  5-  46) 

CALL  PLOTOC !  P , P A ) 

WRITE! 1 , 47) 

WR  I  TE  !  5»  47) 

CALL  PLOTOC !P,  ASNC ) 

WRITE! 1. 48) 

\  WR  I  TE  (  5-  48 ) 

\CALL  PLOTOC  (P,  AOQC  ) 

' WRITE! 1  -  49) 

WRITE!  5.  49) 

CALL  PLOTOC (P, AF I ) 

45  FORMAT! 1H1. 50X,  'OPERATING  CHARACTERISTIC  CURVE') 

46  FORMAT ( 53X <  'PROBABILITY  OF  ACCEPTANCE') 

47  FORMAT! 1H1, 50X,  'AVERAGE  SAMPLE  NUMBER') 

48  FORMAT < 1H1 , 50X,  'AVERAGE  OUTGOING  QUALITY') 

49  FORMAT! 1H1, 50X,  'AVERAGE  FRACTION  INSPECTED') 

GO  TO  898 
556  CONTINUE 

WRITE! 1, 555 ) , LNC 
WR  I  TE !  5.  555),  LNC 

555  FORMAT! 10X,  'LIMIT  NUMBER  FOR  REDUCED  INSPECTION  IS 
898  RETURN 

END 


,  19) 
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OOOl 

1  UNO  1  If  IN  |HN 

P  h 

0002 

PM-Pttl'l 

ooo:i 

G  1.  i' 

ooo* 

Jf  U:N ,  G  i  .  .  i .  - 

•  »  C  '  •  « 

ooo:< 

13  1  XP  '  !‘N  > 

OOO/j 

Z  EXP<  I'N  > 

0007 

IP  (  K.  KG.  0  )  in 

TURN 

ooon 

DU  10  1  IK 

0009 

7- Z  -frPN/ l 

0010 

1  0 

n=--B+z 

001  1 

RETURN 

0012 

1 

B  Q**N 

00 1 13 

Z=Q*  *N 

0014 

IP  (K.  EG.  OIKE 

f  URN 

001  5 

DU  20  I  -  1  <  K 

0016 

Z--Z*  <  N-I  «■  1  ItfPOU*! 

0017 

20 

H-B  +  Z 

0018 

RETURN 

0019 

END 

.SUHROUI  INI  f  NV'  R  I  (NR,  AO,  HQ  ) 


(.)<  *(.'  I 

ooo/ 
o  oo:i 

OO'M 
0005 
ooo<s 
000/ 
0008 
0007 
0010 
OOl  1 
0012 
0018 
OOM 
00  1  0 
00 1  6 
001  7 
0018 
00  IV 
0020 
002  1 
0022 
0023 
008-1 
0020 
0026 
002/ 


r  «  *»  i  HI  6  SUBROUTINE  INVERTS  1  HF  MATRIX  QV  GAUSSIAN  ** 

C  ••  1  I  I  INI  NA  I  IUN.  «v 

< 

REALMS  AQ  <  20,  20  )  ,  BQ<20,  20),  7  I AMP,  A7EMF 
DO  15  1=1, NR 
DO  10  J=  1 ,  NR 
10  8Q(I,J)-0.0 

15  BQ ( I , I ) - 1 . O 

DO  35  1=1, NR 
ZTAMP=AG ( 1,1) 

DO  20  J  - 1 ,  NR 
AO  <  I .  J  )  -  AO  (  I .  J  >  /  Z 1  AMP 
30  13 Q  (  I ,  J  )  -  BO  (  I ,  J  )  /  7  T  AMP 

DU  30  I ]  1 ,  NR 

IF  I  1.1:0.11)  GO  TO  30 
A  TEMP  AO (II,  1  ) 

DO  85  J  1  ,  NR 

AO ( I  I ,  J )  AO ( I  I , J ) -AQ ( I ,  J ) *ATEMP 
25  DO (II,  M  -HO < I  I,  J) -BG< I,  J>*ATEMP 

00  FONT  IN1  >!• 

05  CONI  I  no; 

RF  TURN 
FND 


0(KH 

0002 

0<>ni 
0004 
OOO', 
0006 
OOO  7 
OOOQ 
0009 
0010 
001  1 
0012 
001  3 
OOH 
001  0 
0016 
0017 
00113 
OO  1  9 
0020 
0021 
0022 
0023 
0024 
0020 
0026 
0027 
0020 
0029 
0030 
0031 
0032 
0033 
0034 
0035 
0036 
0037 
0038 
0039 
0040 
0041 
0042 
0043 
0044 


•  I'RRUUl  INI.  (‘LljniCIP  T) 

(. 

(  »*•  IMIS  SUBROUTINE  PI  GTS  OC.  ASM.  AO«  AND  AF 1  CURVES  FOR  «« 
O'  IHf  1  Ul  i  RANG!  UF  INCOMING  FRACTION  DEFECT  IVL.  '  r 

c 

D  I  MENS  ION  U  (11).  TOO).  P(  30),  LINE  (  101  ) 

INTEGER  ASTERK,  BLANK, PI  US 
DATA  ASTERK.  BLANK,  PLUS/  '* 

DO  1  K= 1 .  101 

1  LINE < K ) -BLANK 
XMAX=T< 1 > 

XMIN-T ( 1 ) 

DO  2  I  —  1 , 21 

IF! T< I ) .  I  T.  XM1N  i  XM I N  =T (  I  ) 

IF ( T ( I  )  .  G f .  XMAX )  XMAX-T(I) 

2  CONTINUE 

If  <  XMAX.  L  f .  1  )  XM  I  N  --0 
RANGE “XMAX -XMIN 
RG~-R  ANGF  /  1  O 
U<  )>  XMIN 
DO  3  I  2,  1  1 

3  U(])--U(J  11+-RG 

WR  T  IT  (  J  ,  130)  (U<  I  ) ,  II,  11) 

WRITE  (5,  130)  (U(  I),  I  -1,  II  ) 
no  4  K -  1 ,  101,  10 

4  I.  1  NE  (  K  )  -'PLUS 

WRITE!  1,  100)  (LINE  (K) ,  K--1 »  101  ) 

WRITE  (5,  100)  (LINE  (K),  K=l>  101  > 

DO  6  K=l. 21 
DO  5  1  =  1,  101 

5  LINE! I )=BLANK 

KPR INT= 1 00* (T  <K)— XMIN) /RANGE+1 . 5 
LINE! 1 )=PLUS 
L I NE ( KPR I NT ) = ASTERK 

WRITEd.  140)P(K).  (LINE(I),  1=1,  101) 

WRITE  (  5,  1 40)P  (  K  ) ,  (LINE!  I  >,  1  =  1,  101  ) 

6  CONTINUE 

lOO  FORMAT ( 4  X,  'PERCENT  DEFECTIVE 101 A1 ) 

130  FORMAT! 12X,  1 1 (2X, E3. 3) ) 

140  FORMAT ( 1 7X, F4. 3, 101A1) 

RETURN 

e:nd 
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0001 
ooop 
0000 
0004 
000  5 
0006 
0007 
0008 
0009 
0010 
001  1 
0012 
0013 
0014 
0015 
0016 
0017 
0018 
0019 
0020 
0021 
0022 
0023 
0024 
0025 
0026 
0027 
0028 
0029 
0030 
0031 
0032 
0033 
0034 
0035 
0036 
0037 
0038 
0039 
0040 
0041 
0042 
0043 
0044 
0045 
0046 
0047 
0048 
0049 
0050 
0051 
0052 
0053 
0054 
0055 
0056 
0057 


mjhruut  ini  i  nn<aql  ,  n/,  lnc  ) 

c 

THIS  SUBROUTINE  FINDS  T HI  I  1M) 1  NUMBER  a* 

C**  FOR  REPUCEP  I NSPLCT I UN.  ** 

c 

LNC -222 

IF (NZ. LE. 29. AND. AQL. LE.  15)  LNC  =  0 
IF(NZ. GE. 30. AND. NZ. LE. 49. AND. AQL. LE. 10)  LNC  O 
IFINZ. GE. 50. AND. NZ. LE. 79. AND. AQL. LE. 6. 5)  LNC=0 
I F ( NZ . GE. 80. AND. NZ. LE.  129. AND. AQL. LE. 4.  0)  LNC=0 
IF (NZ. GE.  130. AND. NZ. LE.  199. AND. AQL. LE.  2.  5)  LNC=0 
I F  (  N  Z . GE. 200. AND. NZ. LE. 319. AND. AQL. LE.  1.5)  LNC=0 
I F  (  NZ.  GE.  320.  AND.  NZ.  LE.  499.  AND.  AQL.  LE.  1.  O)  LNC=0 
IF (NZ. GE.  500.  AND. NZ. LE. 799. AND. AQL. LE.  0.  65)  LNC=0 
IF(NZ.  GE.  BOO. AND. NZ. LE.  1249. AND. AQL. LE. 0. 40)  LNC=0 
IF(NZ. GE. 1250. AND. NZ. LE. 1999. AND. AQL . LE. 0. 25)LNC=0 
IF (NZ.  GE.  2000. AND. NZ. LE. 3149. AND. AQL . LE. 0.  15>LNC=0 
IE (NZ. GE. 3150. AND. NZ. LE. 4999. AND. AQL . LE.  0.  10)LNC=0 
IF ( NZ . GE. 5000. AND. NZ. LE. 7999. AND. AQL . LE. 0.  065)LNC=0 
IF ( NZ.  GE.  8000. AND. NZ. LE.  12499. AND. AQL. LE.  O.  040>LNC=0 
I F ( NZ . GE.  12500. AND. NZ. LE.  19799. AND. AQL.  LE.  O.  025)LNC=0 
IF (NZ. GE. 20000. AND. NZ. LE. 31499. AND. AQL  LE. O. 015)LNC=0 
IF  (NZ.  GE.  31500.  AND.  NZ.  LE.  47779.  AND.  AQL..  EG.  0.  010)LNC=0 
IF (LNC. EG. O)  GO  TO  10 
IF(NZ. LE. 29) THEN 

IF  (AQ1  .  EQ.  25 ) LNC--2 
IF  (  AQL.  EQ.  40 ) LNC --4 
IF (AQL. EQ. 65 ) LNC— 8 
IF ( AQL. EQ. 100) LNC— 14 

END  I F 

IF  ( NZ.  GE.  30.  AND.  NZ.  LE.  49 >  THEN 
IF ( AQL. EQ. 15) LNC= 1 
IF ( AQL. EQ.  25 ) LNC=3 
IF (AQL. EQ. 40 ) LNC=7 
IF (AQL. EQ.  65)LNC=13 
IF( AQL. EQ.  100)LNC»22 

END  IF 

IF ( NZ.  GE.  50.  AND.  NZ.  LE.  79)  THEN 
IF ( AQL. EQ. 10 ) LNC=2 
IF (AQL. EQ. 15) LNC=3 
IF (AQL. EQ. 25 ) LNC=7 
IF( AQL. EQ.  40)LNC  =  14 
IF ( AQL. EQ. 65 ) LNC=25 
IF(AQL  .EQ.  100 ) LNC=40 

END  I  F 

IF ( NZ. GE. 80. AND. NZ.  LE.  129) THEN 
IF(AQL . EQ. 6. 5)LNC-2 
IE (AQL. EQ.  1 0  )  LNC -4 
IF  (  AQI  .  EQ.  1  5  )  LNC ~7 
IF( AQL. EQ. 25 ) LNC- 1 4 
IF(AQL. EQ. 40 ) LNC=24 
I F  (  AQI  .  EQ.  65  )  LNC  =  42 
IF (AQL . EQ. 100 )LNC^60 

L  ND  I  F 

IF(NZ. GE. 130. AND. NZ. LE. 199) THEN 


0(  j7 
007)0 
0061 
0062 
0062 
0064 
0065 
0066 
0067 
006G 
0069 
0070 
0071 
0072 
0072 
0074 
0075 
0076 
0077 
00/0 
0079 
0080 
0081 
0082 
0083 
0084 
0085 
0086 
0087 
0088 
0089 
0090 
0091 
0092 
0093 
0094 
0095 
0096 
0097 
0098 
0099 
0100 
0101 
0102 
0103 
0104 
0105 
0106 
0107 
0108 
0109 
0110 
Ol  1 1 
01  12 
0113 
0114 


:  i  '  AQl  .  KG.  4. 0  )  I.  N(  , 

It  (  AGl  .  CO.  6.  3)1.  NO  4 
If  (  AQl  .  LO.  10)  LNC  -  / 

II  <  Alii  .  EQ.  1  5  )LNC-  13 
li-  k  AQL.  EQ.  25  ) LNC  25 
T F  (  AQl  .  EQ.  40)1  NO-4? 

IF(AQL.  EQ.  65  )LNC=72 
IF  (AQl  .  EQ.  1 00  )  LNC“  113 

END  I  f 

IF  (  NZ .  GE.  200.  AND.  NZ.  LE.  319>THF-N 
IF (AQL. EQ. 2. 5>LNC=2 
1  f-  (AQL.  EQ.  4.0) LNC=4 
IF  (AQL.  EQ.  6.  5  )LN08 
IF (AQL. EQ.  10  >  LNC= 1 4 
IF (AQL. EQ. 1 5 ) LNC=22 
IF ( AQL. EQ. 25 ) LNC=40 
IF ( AQL. EQ. 40 ) LNC=60 
IF (AQL. EQ. 65 ) LNC= 1 1 5 
I F ( AQL . EQ. 100)LNC=101 

END  IF 

IF ( NZ. GE. 320. AND. NZ. LE.  4991THLN 
IF  ( AQL.  EQ.  1.  5 ) LNC  =  1 
IF (AQL. EQ. 2. 5>LNC=4 
I F ( AQL. EQ. 4>LNC-8 
IF (AQL. EQ. 6. 5)LNC=14 
I F ( AQL. EQ. 10) LNC ”24 
IF ( AQL. EQ. 1 5 ) LNC-39 
IF (AQL. EQ. 25 ) LNC-68 
IF (AQL. EQ. 65 ) LNC- 1 13 
IF  <  AQL. EQ.  1 00 ) LNC= 1 89 

END  I F 

IF ( NZ. GE. 500. AND. NZ. LE. 799 ) THEN 
I F ( AQL. EQ. 1.0) LNC=2 
IF (AQL. EQ. 1. 5)LNC=3 
IF (AQL. EQ. 2. 5>LNC=7 
IF (AQL. EQ. 4.0)LNC=14 
IF ( AQL. EQ. 6. 5>LNC=25 
IF ( AQL. EQ. 10)LNC=40 
IF (AQL. EQ. 15)LNC=63 
IF  (AQL.  EQ.  25>LNC«=110 
IF ( AQL.  EQ. 40>LNC=181 

END  IF 

IF (NZ.  GE. 800. AND. NZ. LE.  1249) THEN 
IF< AQL. EQ. 0. 65 ) LNC— 2 
I F  <  AQL. EQ.  1 . O ) LNC =4 
IF ( AQL. EQ. 1. 5 ) LNC=7 
IF ( AQL.  EQ.  2.  5)LNC=14 
IF (AQL. EQ. 4. 0)LNC=24 
IF ( AQL. EQ. 6.  5)LNC=42 
IF ( AQL. EQ. 10 )LNC=68 
IF ( AQL. EQ. 15) LNC= 1 05 
IF ( AQL. EQ. 25)LNC=101 

LNDIF 

IF (NZ. GE. 1250. AND. NZ. LE. 1999)THEN 
IF (AQL. EQ. 0. 40 ) LNC-2 
I F  <  AQL. EQ. 0.  65>LNC=4 
I F ( AQL. EQ. 1. 0)LNC^7 
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U  1  1  ■  # 
01  »/. 
011/ 
0111) 
011  V 
0 1  20 
0121 
0122 
0123 
0124 
0125 
0126 
0127 
0128 
0129 
0130 
0131 
0 1  32 
0133 
0134 
0135 
0136 
0137 
0130 
0137 
0140 
0141 
0142 
0143 
0144 
0145 
0146 
0147 
0148 
0149 
0150 
0151 
0152 
0153 
0154 
0155 
0156 
0157 
0158 
0159 
0160 
0161 
0162 
0163 
0164 
0165 
0166 
0167 
0168 
0169 
0170 
0171 


1  I  (  AUL.  EQ.  1  .  5  )LNC  1  1 
1  I  (AOL  .  E  Q.  2.  5  >  LNC  -24 
II-  (AOL.  EQ.  4.  OJLNC  -40 
IE  (AOL..  EQ.  6.  5  )  LNC-69 
If-  (  AOL .  EQ.  1  0  )  LNC-  1  10 
IE (AOL. EQ. 15JLNC- 169 
END  If 

I F ( N7 .  GE. 2000. AND. N2. LE.  3149)  TUFN 
IF (AOL. EQ. 0. 25 ) LNC-2 
I F  (  AQL.  EQ.  0.  401LNO4 
I F  <  AQL. EQ. 0. 65 ) LNC-8 
IF (AQL. EQ. 1. 0 ) LNC- 1 4 
I F  <  AQL.  EQ.  1.5)  LNC-22 
IF  (AQL.  EQ.  2.  5)LNO40 
IF (AQL. EQ. 4. 0)LNC-68 
IF (AQL. EQ. 6. 5>LNC=105 
IE (AQL. EQ. 10)LNC=101 


END1J- 
I  F  (  N  Z  . 


ENDIE 
IF  (N7 


GE. 3150. AND. NZ. LE. 4999) THEM 
1 F ( AQL. EQ. 0. 15) LNC- 1 
I F ( AQL. EQ. 0. 25 ) LNC=4 
IF (AQL. EQ. 0. 40 ) LNC=8 
I F ( AQL. EQ. O. 65>LNC=14 
I F ( AQL. EQ. 1.0) LNC- 24 
I F  ( AQL.  EQ.  1. 5)LNC'-38 
IF  (AQL.  EQ.  2.  5)LN0  67 
IF (AQL. EQ. 4.0) LNC" 1 1 1 
I F ( AQL. EQ. 6. 5)LNC= 186 


[ F (N7 .  GE. 5000. AND. NZ. LE. 7999) THEN 
IF (AQL. EQ.  O.  10  1LNC-2 
IF ( AQL. EQ. O. 1 5) LNC =3 
IF ( AQL. EQ. 0. 25)LNC=7 
IF ( AQL. EQ. 0. 40 ) LNC=1 4 
IF<  AQL. EQ. O. 65 ) LNC— 25 
I F ( AQL. EQ. 1 . O ) LNC=40 
IF (AQL. EQ. 1. 5)LNC=63 
IF (AQL. EQ. 2. 5 ) LNC=I 10 
IF (AQL. EQ. 4. 0)LNC=181 

ENDIF 

IF(NZ.  GE. 8000. AND. NZ.  LE.  12499) THEN 
IF( AQL. EQ. 0. 065 ) LNC- 2 
IF ( AQL. EQ. O. 1 O ) LNC=4 
IF ( AQL. EQ. 0. 15) LNC=7 
IF (AQL. EQ. O. 25>LNC=14 
IF (AQL. EQ. 0. 40)LNC=24 
IF ( AQL. EQ. 0. 65>LNC=42 
I F  ( AQL.  EQ.  1. 0>LNO68 
IF (AQL. EQ. 1. 5 ) LNC~ 1 05 
IF (AQL. EQ. 2. 5)LNC= 181 

ENDIF 

IF (NZ. GE. 12500. AND. NZ. LE. 19999)THEN 
I F ( AQL. EQ. 0. 040 ) LNC=2 
IF ( AQL. EQ. 0. 065 ) LNC=4 
IF(AQL. EQ. 0. 1 O ) LNC=7 
I F  (  AQL.  EQ.  0.  1  5 )  LNC:_  1 3 
IF ( AQL. EQ. 0. 25>LNC=24 


i 


O'  ';t 
O  !  /4 

o  i  /:• 

01  70 
oi  // 
oi  /n 
017V 
0180 
0181 
0182 
0183 
0184 
0185 
0186 
0187 
0188 
018V 
0190 
0191 
0192 
0193 
0194 
0195 
0196 
019/ 
0  1  93 
0199 
0200 
0201 
0202 
0203 
0204 
0205 
0206 
0207 
0203 
0209 
0210 
021 1 
0212 
0213 


END  11 
1 1  •'  <  N  7  , 


END  I  V 
IF (NZ. GE. 

IE  - 

ir 

ir 

IF 

IF 

IF 

IF 

IF 

IF 

FNDIF 
IF (NZ. GE 
IF 
IF 
IF 
IF 
IF 
IF 
IF 
IF 
IF 
IF 

ENDIF 

CONTINUE 

RETURN 

END 


*  AGI.  .  EG.  O.  40)1.  N(  40 
'  AUL  .  EG.  0.  65  /  LNC  69 
(AQL.  EG.  1.  OH.NC-  110 
:  (  AQL.  EG.  1 . 5)1. NC  169 

.  20000.  AND.  NZ.  i_E.  3149V 
'■  (  AQL.  EG.  0.  025) LNC ~2 
•  (AQL.  EQ.  0.  040)1  NC=-4 
-(AGL. EG. 0. 065 ) LNC=8 

-  ( AQL. EG. O.  10) LNC- 14 

-  (AQL. EG. 0.  15)LNC=22 
-(AGL. EG. 0. 25 ) LNC=40 

-  ( AGL. EG. O. 40 ) LNC=68 

:  ( AQL. EG. 0. 65 ) LNC= 115 
'  (AGL. EG.  1. 0)LNC=1B1 


.  3 1 500. AND. NZ . LE.  49999 ) THEN 
(AQL. EG. 0.  015) LNC  =  1 
(AQL. EQ. O. 025)LNC=4 
(AQL.  EG.  0.  040) l-NC ~B 
(AQL. EG. 0. 065 ) LNC= 14 
(AQL. EG. 0. 10)LNC-24 
(AGL.  EQ.  O.  15) *-NC=38 
(AQL. EG. 0. 25 ) LNC=67 
( AQL. EG. O. 40 ) LNC=1 1 1 
(AQL. EG. O. 65) LNC- 186 


.  50000)' 
(AQL. EQ 
(AQL. EG 
(AQL. EG 
(AQL. EG 
(AGL. EQ 
(AQL. EG 
(AQL. EG 
(AQL. EG 
(AQL. EQ 
(AQL. EQ 


THEN 

. O. O 1 O ) LNC=2 
. O. 015)LNC=3 
. O. 025 )LNC=7 
. O. 040 ) LNC= 1 4 
. 0. 065 ) LNC=25 
. 0. 10)LNC=40 
. O. 1 5 ) LNC=63 
.0.  25)LNC=HO 
. O. 40)LNC=181 
. 0. 65)LNC=301 


*v*  •J'  V-  V* 


v! 


Program  PR0BS1.F0R 


Quo  i 

0002 

UUO.< 

000-1 

ooo:> 

000,'. 

0007 

000(1 

0009 

OO  V  o 
001  1 
0012 
OO  1  3 
001-1 
OO  1  0 
00  1  A 
OOl  / 
OOl  B 
0019 
0020 
0021 
oop;' 

0023 
002-1 
0025 
0026 
002/ 
0028 
0029 
0030 
0031 
0032 
0033 
003-4 
0035 
0036 
0037 
0038 
0039 
0040 
0041 
0042 
0043 
0044 
004  5 
0046 
004  7 
0048 
0049 
0050 
0051 
0052 
0053 
0054 
0055 
0056 
0057 


SUMR0U1  INh  PROBSllNN,  f-  •  .  U/  i  Li 

r  o  v  f.  i.  *•*«■<<  *********#**********■"*** 

i  n  i  r>  SUBROUT  I  NT  COMPUTT  ruMi.n  ATTVf  RTMUWAI. 

(  PPOBABIL  IT1ES  «... 

IN  I  EGER  C 

DOUBLE  PRECISION  SUMLOG 

( . 

uOMMQN/BLK7 /SUMLOG  <  4000 ) 

COMMDN/BL K8/N 

O  ■  - l .  -v- 

C  BINOMIAL  PROB.  WHEN  C=0 

C  *  o  *  ****  t!  ***«**♦<(* 

( SUMS  O* *NN 

TT  (C.  [tq.  o)  GOTO  -15 

***»***»♦********»»«'***»»*»«»«'*»**♦♦»*♦«»*«  K  ** 

C  AVOID  RECOMPUTING  SUMLOGdl'S  ALREADY  IN  MEMORY 

******** 

IP  (N  NN)  10,  25,  25 
10  N  M-*  1 

C  COMPUTE  N  SUMLQGS-EQUIVALENT  TO  N -FAC TOR  I AL 

C******* 

IF  (M. GT.l)  GOTO  15 
SUMLOG ( 1 1^0. 

IF (NN. LE.  1  )  GOTO  25 

M-2 

15  DO  20  I=M, NN 

SUMLOG ( I ) =DLOG 1 O ( DFLOAT ( I ) ) +SUMLOG ( I - I ) 

20  CONTINUE 

C ****************************************************** 

C  COMPUTE  C  SUMS-EQU I VALENT  TO  SSUM  OF  PROB.COMPIN. 

C  I.E.  CUMULATIVE  BINOMIAL  DISTRIBUTION  COMPUTATION 

C ****************************************************** 

25  IF(NN.  GT.N)  N=NN 

C ****************************************************** 

C  DETERMINE  BEST  NUMBER  HANDLING  LOOP 

C ************************************************ ****** 

IF  (NN.GT.300)  GOTO  35 
DO  30  K=l.  C 

CSUMS= 10.  **( SUMLOG < NN ) -SUMLOG < NN-K ) -SUMLOG (K) ) 

1  *P**K*Q**  (  NN-K  )  +-CSUMS 

30  CONTINUE 
GOTO  45 

C  ****************************************************** 

C  LOOP  FOR  LARGE  EXPONENTS 

C  ****************************************************** 

35  DO  40  K= 1 ,  C 

CSUMS= 1 0. ** ( SUMLOG ( NN ) -SUMLOG ( NN-K ) -SUMLOG ( K  > 

1  +K*DL  OG 1 0  <  DBLE ( P  )>  +  < NN-K ) *DLOG 1 0  <  DBLE  <  Q ) ) ) +CSUMS 

40  CONTINUE 
C 

45  B XI  EC  -  CSUMS 
RETURN 
END 


Program  PR0BD1.F0R 
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